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ABSTRACT 


Autonomous  object  recognition  is  an  active  area  of  interest  for  military  and 
commercial  applications:  Given  an  input  image  from  an  infrared  or  range  sensor,  find 
interesting  objects  in  those  images  and  then  classify  those  objects.  In  this  work,  automatic 
target  recognition  of  ship  types  in  an  infrared  image  is  explored.  The  first  phase  segments 
the  original  infrared  image  in  order  to  obtain  the  ship  silhouette.  The  second  phase 
calculates  moment  functions  of  those  silhouettes  that  guarantee  invariance  with  respect  to 
translation,  rotation  and  scale.  The  third  phase  applies  those  invariant  features  to  a 
backpropagation  neural  network  and  classifies  the  ship  as  one  of  the  five  types.  The 
algorithm  was  implemented  and  experimentally  validated  using  both  simulated  three- 
dimensional  ship  model  images  and  real  images  derived  from  video  of  an  AN/AAS-44V 
Forward  Looking  Infrared  (FLIR)  sensor. 
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The  algorithms  and  computer  programs  developed  in  this  research  were  not 
exercised  for  all  possible  cases  of  interest.  While  every  effort  has  been  made,  within  the 
time  available,  to  ensure  that  the  programs  are  free  of  computational  and  logic  errors, 
they  cannot  be  considered  validated.  Any  application  of  these  programs  without 
additional  verification  is  at  the  risk  of  the  user. 
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1.  INTRODUCTION 


A.  AUTOMATIC  TARGET  RECOGNITION 

Automatic  target  recognition  (ATR)  is  a  technological  discipline  that  deals  with 
the  understanding,  design,  development,  and  production  of  techniques  and  hardware  for 
the  classification  of  objects  of  interest  as  they  are  sensed  by  remote  means,  either  actively 
or  passively.  During  the  past  several  decades,  numerous  attempts  have  been  made  to 
create  such  systems  [Ref  1-6],  However,  progress  in  ATR  has  been  slow  [Ref  7] 
because  some  new  problems  have  appeared.  For  example,  the  vision  problem  pushes  the 
fields  of  artificial  intelligence,  neural  networks,  microelectronics,  sensors,  and  computer 
science  to  their  limits. 

In  any  pattern  recognition  application,  it  is  important  to  select  features  that 
adequately  and  uniquely  describe  the  objects  to  be  recognized.  Moreover,  the  features 
associated  with  an  object  should  be  invariant  with  respect  to  the  position,  rotation,  and 
scale  of  that  object  in  the  field  of  view.  Thus  the  ideal  recognition  system  is  robust  to 
orientation  variations,  scale  variations  and  boundary  perturbations  [Ref  8]. 

Our  proposed  approach  is  to  use  the  moment  invariants  [Ref  9]  for  the  set  of 
features  to  quantify  the  object.  The  thesis  reports  the  mathematical  foundation  of  two- 
dimensional  moment  invariants  and  shows  that  recognition  schemes  based  on  them  could 
be  truly  position,  size  and  orientation-independent.  Since  the  moments  are  global 
features,  application  of  such  a  feature  space  is  limited  to  images  with  minimal 
background  and  scenes  containing  only  one  object.  Ships  on  the  open  sea  are  appropriate 
for  such  feature  spaces  [Ref  10].  The  moment  invariants  are  used  to  construct  a  feature 
vector  of  low  dimension,  and  recognition  is  performed  using  this  feature  vector  applied  to 
a  trained  artificial  neural  network  classifier. 
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B.  APPLICATIONS  OF  AUTOMATIC  OBJECT  RECOGNITION 


Automatic  object  recognition  has  diverse  applications  in  numerous  fields  of 
science  and  technology  and  is  permeating  many  aspects  of  military  and  civilian 
industries.  It  is  popular  within  the  field  of  robotic  vision  because  of  the  limited  domains 
and  the  controllability  of  the  environment  in  which  they  are  used  [Ref.  7]. 

In  military  applications  and  specifically  naval  applications,  electro-optic  and 
infrared  sensors  have  been  connected  to  weapon  systems.  In  several  cases  an  ATR 
algorithm  is  responsible  for  discriminating  a  target  from  a  non-target  object,  enabling  the 
possible  target  destruction.  Other  systems  address  classification  tasks  where  targets  types 
are  determined.  Another  example  is  a  Forward  Looking  Infrared  (FLIR)  sensor  combined 
with  image-processing  hardware  for  discriminating  tanks  from  trucks,  bushes,  and  other 
environmental  objects. 


C.  PROJECT  GOALS 

In  this  thesis,  a  fast  and  robust  system  is  presented  that  classifies  ships  seen  from 
an  arbitrary  viewpoint  and  range  in  three-dimensional  space.  The  approach  is  concerned 
with  segmented  rigid  bodies  viewed  without  occlusion  from  other  objects.  However,  self¬ 
occlusion  due  to  change  of  viewpoint  is  allowed. 

Although  object  separation  from  background  is  a  challenging  task  in  general,  our 
application  can  be  carried  out  with  relative  ease.  This  is  the  case  because  a  ship  has 
usually  a  clear  contrast  with  the  background  in  Forward  Looking  Infrared  (FLIR) 
imagery.  This  greatly  simplifies  ship  classification. 

Our  approach  is  model-based,  meaning  that  the  kinds  of  objects  to  be  recognized 
are  known  in  advance  and  can  be  summarized  in  a  set  of  models.  The  specific  model 
database  we  have  implemented  contains  five  classes  of  ships:  destroyer,  frigate,  aircraft 
carrier,  research  ship  and  merchant  ship.  This  database  was  used  in  the  training  phase. 
For  each  ship  model  and  viewpoint,  a  silhouette  was  extracted  and  a  moment-invariant 
signature  calculated  consisting  of  a  twelve-element  feahire  vector. 
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Then  for  a  ship  image  of  unknown  type,  we  compute  its  signature.  Classification 
is  done  using  an  artificial  neural  network.  The  neural-net  classifier’s  generalization 
capabilities  are  used  to  group  the  moment-invariant  signatures,  corresponding  to  different 
views  of  an  object,  into  a  single  ship  type  class. 

The  proposed  scheme  is  summarized  in  Figure  1  below. 


Figure  1 :  Processing  scheme  of  moment  invariant  recognition. 


This  thesis  has  been  organized  as  follows:  in  Chapter  II  we  present  techniques  for 
feature  extraction  followed  by  mathematical  foundations  of  moment  invariants.  In 
Chapter  HI  we  present  an  artificial  neural  network  as  a  classifier  and  the  backpropagation 
learning  rule.  Chapter  IV  details  the  input  image  database  and  the  training  and  testing  of 
our  system.  Chapter  V  summarizes  the  results  from  the  experimentation  using  simulated 
images  from  three-dimensional  ship  models  and  real  ship  images  from  FT, TP  sensors. 
Chapter  VI  contains  concluding  remarks. 
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11.  FEATURE  SELECTION  AND  MOMENT  INVARIANTS 


A.  OVERVIEW 

An  old  adage  says  “Good  Features  make  Good  Recognizers”  [Ref.  11].  This  is 
true  whether  your  recognizer  is  using  an  artificial  network  or  a  statistical  based  decision 
mechanism.  So  our  project  paid  careful  attention  to  feature  extraction  from  ship  images. 

Many  approaches  have  been  advocated  for  features  in  automatic  target 
recognition.  Present  methods  can  be  categorized  as  either  global  or  local.  Global  methods 
use  global  features  of  an  object  boundary  or  of  an  equivalent  representation.  Such 
techniques  are  the  Fourier  descriptors  (FD)  [Refs.  1,  5],  moments  [Ref  12]  and 
autoregressive  models  [Ref  13].  Local  methods  use  features  such  as  critical  points  [Ref 
14]  or  high-resolution  pursuit  (HRP)  [Refs.  8, 15]. 

For  global-based  approaches,  there  is  a  wide  variety  of  published  literature  similar 
to  the  approach  described  herein.  Global  methods  have  the  disadvantage  that  a  small 
distortion  in  a  section  of  a  boundary  of  an  object  will  result  in  changes  to  all  global 
features. 

One  early  work  is  Dudani  et  al  [Ref  16],  which  used  moment  invariants  for 
feature  extraction  and  a  probabilistic  approach  for  the  classification  of  airplanes.  Dudani 
used  six  different  aircraft  types  and  the  images  were  based  on  physical  models.  His 
training  set  was  based  on  over  3000  images  taken  in  a  140°  by  90°  sector.  The  testing  set 
contained  132  images  (22  images  of  each  of  the  six  classes)  obtained  at  random  viewing 
aspects.  The  classification  accuracy  achieved  in  this  six-class  problem  was  95%. 

Later,  Wallace  and  Wintz  [Ref.  17]  propose  a  technique  similar  to  Dudani’s  with 
Fourier  Descriptor  (FD)  of  the  silhouette  boundary  as  features.  The  Fourier  descriptor  is 
one  method  of  describing  the  shape  of  a  closed  figure.  Wallace  and  Wintz  used  a  graphics 
program  to  test  their  algorithm  implementing  three-dimensional  models  for  six  different 
aircrafts.  The  graphics  program  approximated  each  airplane  by  using  50-100  planes. 
Again,  the  evaluation  was  done  using  a  randomly  selected  set  and  comparing  to  the 
library  of  projections.  However,  they  used  only  143  projections  for  training  (9.9  times 
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less  than  that  used  by  Dudani  et  al)  and  the  aircraft  outlines  were  taken  from  a  sector  of 
180°  by  180°.  Wallace  and  Wintz  considered  a  bigger  sector  trying  to  avoid  Dudani’s 
approach,  because  if  we  delete  the  angles  near  the  front  view  and  rear  view  of  the  aircraft 
the  problem  is  much  easier:  shapes  vary  much  more  with  slight  rotations  when  viewed 
almost  edgewise  [Ref.  17].  The  maximimi  classification  accuracy  achieved  by  Wallace 
and  Wintz  was  88.0%. 

Reeves  et  al  [Ref.  18]  presents  a  geometrical-moment  approach  using  moments  of 
the  image  that  are  normalized  with  respect  to  scale,  translation  and  rotation.  They  call 
them  “standard  moments”.  The  experiments  described  there  were  based  on  the  same 
software  used  by  Wallace  and  Wintz.  They  also  used  the  same  six  types  of  airplanes,  the 
same  training  set  and  the  same  testing  set.  However,  they  have  chosen  the  moment 
feature  representation  because  Fourier  descriptors  (FD)  are  particularly  sensitive  to 
perturbations  in  the  object  boundary.  For  example,  the  FD’s  for  the  image  of  a  disk  differ 
greatly  from  those  for  a  disk  with  a  tiny  wedge  missing.  Reeves  et  al  used  two 
classification  criteria:  the  minimum  Euclidean  distance  and  the  minimum  Euclidean 
distance  after  “variance  balancing”.  The  best  classification  result  was  93%. 

More  recent  work  of  Khotanzad  [Ref.  19]  used  global  features  derived  from 
complex  orthogonal  Pseudo-Zemike  Moments  (PZM).  Khotanzad  tested  the  performance 
of  PZM  by  recognizing  26  uppercase  English  characters  (A  to  Z),  typed  and  handwritten. 
The  database  contained  624  images  corresponding  to  24  images  per  character.  These 
images  were  generated  with  arbitrarily  varying  scales,  orientations,  and  translations.  The 
available  samples  were  divided  into  halves.  The  first  half  was  used  for  training  and  the 
second  for  testing.  There  were  12  training  images  and  12  testing  images  per  character. 
His  neural  network  classifier  formed  by  45  input  nodes,  26  output  nodes  and  40  hidden 
nodes  got  100%  of  classification  accuracy. 

Systems  using  local  features  perform  well  in  the  presence  of  noise,  distortion  or 
partial  occlusion.  The  effects  on  an  isolated  region  of  the  contour  alter  only  the  local 
features  associated  with  that  region,  leaving  all  the  other  local  features  unaffected. 
However,  the  choice  of  representative  local  features  is  not  trivial  and  the  recognition 
process  based  on  local  features  is  more  computationally  intensive  and  time  consuming 
[Ref  2]. 
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B.  MOMENT  INVARIANTS 


Moment  invariants  are  a  reliable  and  versatile  way  to  construct  a  feature  vector  of 
low  dimension  as  the  basis  for  the  neural-network  classifier.  Moments  have  been  used  as 
pattern  features  in  a  number  of  applications  [Ref  9,  20]  to  recognize  two-dimensional 
image  patterns. 

The  regular  moments  mpq  of  a  digital  image  pattern  represented  by  f(x,y)  are 
defined  as: 


p,g=ox2,...  a) 

jc  y 

Hu  [Ref  9]  first  introduced  moments  as  image-recognition  features.  Using 
nonlinear  combinations  of  normalized  central  moments,  he  derived  seven  invariant 
moments,  which  have  the  desirable  property  of  being  invariant  under  image  translation, 
scaling  and  rotation.  The  classic  central  moments  that  have  the  property  of  translation 
invariance  are: 


p,g=o.i,2.. 


(2) 


^^0  _  ,  -  ^01 


where  x  = -  and  y  =  - 


m 


00 


m 


00 


Hu  discovered  these  moments  Ml,  M2,  ...,  M7,  are  invariant  under  translation 
and  rotation: 


Ml  =  )i20+  M02 

(3) 

M2  =  (p20  -  lk>2)^  +  4)0,^  1 1 

(4) 

M3  =  (P30-  3pi2)^  +  (3p2i  - 

(5) 

M4  =  (p3o+  M-u)^  +  (3p2i  +  1^03)^ 

(6) 
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M5  =  (|i30-  3^12X^30  +  I^i2)[  (^30+  -  3(|J2i  +  Hos)^  ] 

+  (3|a2l  -  )X03X  fX21  +  |^03)[  3(^30  +  1^12)^  “  (^21  +  M03)^  ]  (7) 

M6  =  (^20  -  M'02)[  ([^30  +  M-12)^  -  (^21  +  ^03)^  ]  +  4^i  i(|J.30  +  Hl2)(  ^21  +  R03)  (8) 

M7  =  (3)121  -  M03)(M-30+  M.12)[  (^30+  |J-12)^  “  3(  II2I  +  M03)^  ] 

-  (1^30-  3)Xi2)  (1^21  -  ^03)[  3(^30  +  fil2)^  -  (1^21  +  1^03)^  ]  (9) 

The  functions  Ml  through  M6  are  invariant  under  rotation,  reflection,  or  a 
combination  of  rotation  and  reflection.  This  property  helps  to  simplify  the  range  of  all 
distinct  views  of  the  ships  as  explained  in  section  D  of  chapter  IV. 

The  above  moments  can  be  normalized  to  become  invariant  under  a  scale  change 
by  using  the  radius  of  gyration  r  of  a  planar  pattern  [Ref.  16]: 

r  =  (p2o+Po2)'^^  (10) 

The  radius  of  gyration  for  a  particular  object  from  a  particular  angle  of  view  is 
directly  proportional  to  the  size  of  the  image  or  inversely  proportional  to  the  distance  B 
of  the  object  along  the  optical  axis: 

(P2o  +  Po2)^^^  B  =  constant  (1 1) 

Therefore,  the  radius  of  gyration  r  can  normalize  the  moment  functions  M2 
through  M7  to  obtain  size  invariance,  what  Hu  called  the  “normalized  central  moments”: 


Ml  ’  —  (p2o  |-k)2)'^^  B  —  r  B 

(12) 

M2’=M2/r^ 

(13) 

M3’  =  M3/r^ 

(14) 

M4’  =  M4  /  r® 

(15) 

M5’=M5/r’^ 

(16) 

M6’=M6/r® 

(17) 
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(18) 


M7’  =  M7/r’^ 

C.  FEATURE  VECTOR 

The  above  moments  of  an  object  can  be  computed  for  both  the  image  boundary 
and  the  solid  silhouette.  Minute  details  such  as  the  shape  of  the  stacks  of  a  ship  are  better 
characterized  by  the  moments  from  the  boundary.  Gross  structural  features  of  the  ship  are 
better  characterized  by  moments  derived  of  silhouette;  also,  these  moments  are  less 
susceptible  to  noise  [Ref  16]. 

In  our  system,  two  sets  of  six  moment  invariant  functions  (M2’,  M3’,  M4’,  MS’, 
M6’  and  M7’),  six  from  the  boundary  and  six  from  the  silhouette,  were  computed.  As  the 
distance  B  of  the  object  along  the  optical  axis  was  not  known,  the  Ml’  component  was 
not  used.  The  twelve-component  feature  vector  was  sent  to  the  neural  network  classifier 
for  the  recognition  phase. 
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III.  THE  ARTIFICIAL  NEURAL  NETWORK  (ANN)  CLASSIFIER 


A.  MOTIVATION 

Over  the  past  few  years,  an  explosion  of  interest  in  ANN  models  and  their 
applications  has  occurred  [Ref  21,  22,  23].  ANNs  posses  a  number  of  properties  which 
make  them  particularly  suited  to  complex  classification  problems  [Ref  22,  25,  26]. 
Unlike  traditional  classifiers,  ANN  models  can  examine  numerous  competing  hypotheses 
simultaneously  using  massive  interconnections  among  many  simple  processing  elements. 
In  addition,  ANNs  perform  extremely  well  under  noise  and  distortion. 

The  implementation  of  a  model-based  target  recognition  scheme  using  ANNs 
seems  to  be  attractive.  First  of  all,  ANNs  provide  their  own  way  to  represent  the 
knowledge  that  they  store  [Ref.  27].  In  addition,  the  complexity  and  the  computational 
burden  increase  slowly  as  the  number  of  data  models  increases. 

Although  ANN’s  performance  is  excellent,  many  researchers  still  criticize  ANNs 
because  they  can  require  much  training  time  before  they  can  perform  a  specific  task. 
However,  in  our  automatic  target  recognition  classifier,  the  recognition  phase  is  of  far 
more  importance  and  it  must  run  as  quickly  and  accurately  as  possible;  the  training  phase 
can  be  performed  off-line. 

In  this  thesis  specifically,  a  three-layer  perceptron  neural  network  [Ref  28] 
trained  with  the  backpropagation  learning  rule  [Ref  29]  was  implemented.  In  this 
scheme,  expensive  storage  of  a  multiview  database  is  not  needed  since  during  training  the 
neural  net  extracts  all  the  relevant  information  from  the  library.  Also,  due  to  the 
generalization  capability  of  the  neural  net,  good  results  can  be  obtained  even  with  a  small 
number  of  views  in  the  library. 
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B.  MULTILAYER  ANN  AND  THE  BACKPROPAGATION  RULE 


Artificial  neural  networks  were  developed  by  modeling  a  biological  neuron.  A 
generic  neuron  is  formed  by  “cell  body”,  “dendrites”  and  “axon”.  The  electrical  signals 
arrive  in  a  neuron  by  the  dendrites  and  are  passed  to  the  cell  body  where  they  are  added. 
If  a  threshold  is  achieved,  the  neuron  is  activated  and  the  information  is  passed  to  the 
axon.  The  axon  is  the  transmission  line  of  the  neuron.  The  axon  will  pass  the  information 
to  the  chemical  synapses  connections.  The  learning  process  will  be  responsible  for 
increasing  the  synaptic  strength,  which  measures  the  degree  of  coupling  between  two 
neurons. 

Many  neuron  models  appear  in  the  literature.  The  beginning  of  the  development 
of  neural  network  models  is  related  to  the  paper  of  Warren  McCulloch  and  Walter  Pitts 
published  in  1943.  They  studied  the  implementation  of  logical  functions  using  artificial 
neurons.  The  mathematical  model  of  the  neuron  proposed  by  McCulloch-Pitts,  shown  in 
Figure  2  below,  assumes  the  function  realized  by  the  cell  body  as  being  a  “step  function” 
applied  to  the  summation  of  the  weighted  inputs.  The  weights  control  the  importance  of 
each  input. 


Neuron 


Figure  2:  McCulloch-Pitts  Model. 


The  McCulloch-Pitts  model  applied  to  a  single  layer  of  neurons  (perceptrons) 
cannot  solve  problems  where  the  inputs  cannot  be  linearly  separated.  The  PDP  group  in 
their  collection  of  papers  [Ref  28]  proposed  modifications  to  the  previous  model.  The 
step  function  was  replaced  by  a  function  that  is  monotonic,  differentiable  and  smooth 
(often  implemented  by  a  sigmoid).  The  learning  algorithm  used  is  “backpropagation”. 
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Artificial  neural  networks  (ANN)  are  specified  by  the  topology  of  the  network, 
the  characteristics  of  the  nodes  (neurons)  and  the  learning  algorithm.  The  topology  of  a 
multilayer  ANN  is  a  structured  hierarchical  layered  network  as  shown  in  Figure  3  below: 


Figure  3:  Multilayer  neural-network  graph. 

It  consists  of  several  layers  of  nodes,  and  usually  an  input  layer  and  an  output 
layer.  Between  the  input  layer  and  the  output  layer,  we  have  one  or  more  “hidden”  layers 
of  nodes.  Hidden  nodes  often  represent  domain  knowledge  useful  for  solving  recognition 
tasks  [Ref.  27].  Generally,  each  node  in  one  layer  is  interconnected  with  all  the  nodes  in 
adjacent  layers  with  connections  (synapses).  Each  connection  is  associated  with  a  weight, 
which  measures  the  degree  of  interaction  between  the  corresponding  nodes. 

A  general  L-layered  feed- forward  artificial  neural  network  consists  of  No  input 
nodes  and  Nl  output  nodes.  The  number  of  nodes  in  the  hidden  layers  is  Nk  for  l<k<L-l. 
In  this  notation,  the  input  layer  is  not  counted  as  a  layer.  So  an  L-layer  feed-forward 
artificial  neural  network  has  L-1  hidden  layers  and  the  Lth  layer  is  the  output  layer.  In 
this  thesis,  we  implemented  a  2-layered  (L=2)  feed-forward  artificial  neural  network  with 
No=12  (the  moment  invariant  feature  vector)  and  N2=5  (five  ship  types).  The  number  of 
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hidden  nodes  was  found  in  order  to  maximize  the  neural  net  performance,  as  detailed  in 
Chapter  V. 

The  algorithms  for  multilayer  ANN  processing  can  be  divided  into  two  phases: 
retrieving  and  learning.  In  the  retrieving  phase  of  the  algorithm,  information  flows  from 
the  input  layer  through  the  hidden  layers  to  the  output  layer.  The  nodes  update  their  own 
activation  values  based  on  the  system  dynamics.  In  the  learning  phase,  modification  of 
the  weights  corresponding  to  the  connection  edges  takes  place.  In  this  thesis,  the  popular 
backpropagation  rule  [Ref  28]  learning  algorithm  is  used.  This  algorithm  performs 
supervised  learning;  in  each  step  it  adjusts  the  connection  weights,  minimizing  the  mean- 
square  error  between  the  target  value  (the  desired)  and  the  output  value  (the  actual)  if  the 
network. 

During  the  retrieving  phase,  we  present  continuous  valued  input  data  xi,  X2,  ..., 
Xno  called  exemplar  patterns  and  the  corresponding  desired  output  data  ti,  tz, . . .,  tnL  called 
target  patterns.  Input  data  are  propagated  forward  through  the  network,  which  computes 
the  activation  value  for  each  node,  until  the  output  layer  is  reached. 

The  learning  phase  involves  a  backward  pass  through  the  network  during  which 
the  error  signals  produced  at  the  output  layer  are  passed  to  each  node  in  the  network  and 
appropriate  weight  changes  are  made.  For  each  weight,  the  gradient  of  the  output  error 
with  respect  to  that  weight  is  computed.  The  weight  is  changed  in  the  direction  that 
reduces  the  error. 
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IV.  EXPERIMENT  DESCRIPTION 


A.  PROGRAMMING  ENVIRONMENT 

Our  Automatic  Recognition  Algorithm  was  based  on  programs  written  in 
MATLAB  5.3.0  from  MathWorks.  MATLAB  is  a  complete  computing  environment  for 
the  interactive  analysis  and  visualization  of  data,  integrating  an  array-oriented  language 
with  mathematical  analysis  and  graphical  display  techniques.  The  neural-network 
programs  used  in  this  thesis  were  implemented  using  functions  from  the  MATLAB 
Neural  Network  Toolbox.  The  three-dimensional  model  and  all  image  analysis  were 
performed  using  the  MATLAB  Image  Processing  Toolbox. 

Although  MATLAB  is  an  interpretative  language,  it  is  possible  to  translate  all 
MATLAB  source  codes  into  C  code  and  create  executable  files  using  the  MATLAB  C 
Compiler.  Consequently,  our  Target  Recognition  System  could  be  used  in  a  real-time 
application. 

B.  THE  THREE-DIMENSIONAL  SHIP  MODEL  DATABASE 


1.  The  Three-Dimensional  Ship  Modeling 

This  section  will  describe  our  implementation  of  three-dimensional  ship  models. 
The  three-dimensional  wire-frame  models  represent  a  graphics  object  by  connected 
polygons  or  faces.  The  model  is  defined  by  specifying  the  coordinates  of  the  vertices  of 
each  polygon  and  then  specifying  the  faces  by  connecting  the  specified  vertices  in  a 
specific  order.  This  three-dimensional  modeling  was  based  on  a  MATLAB  fimction 
called  “patch”  (see  “findInputSet.m”  in  the  Appendix  A). 

Five  ship  types  were  chosen  to  be  included  in  the  recognition  class  and  therefore 
be  modeled:  namely,  an  aircraft  carrier,  a  fngate,  a  destroyer,  a  research  ship  (Point  Sur), 
and  a  merchant  ship.  With  these  five  types,  it  was  possible  to  address  a  typical  scenario  at 
sea,  where  we  can  find  military  ships,  small  civilian  ships  and  big  merchant  ships. 
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The  three-dimensional  wireframe  model  for  the  aircraft  carrier  was  based  on  a 
1:1800  scaled  drawing  of  the  Carl  Vinson  aircraft  carrier  (Nimitz  Class)  [Ref  30].  This 
drawing  is  shown  in  Figure  4  and  a  picture  of  the  Carl  Vinson  aircraft  carrier  is  shown  in 
Figure  5.  The  model  was  implemented  manually  since  no  CAD  model  was  available.  It 
was  formed  of  45  vertices  and  34  planes  (see  the  “aircraft.m”  program  in  Appendix  A). 
Figure  6  shows  the  model  from  four  view  angles. 


Figure  4:  Scaled  drawing  of  the  Carl  Vinson  aircraft  carrier  [From  Ref.30] 


Figure  5:  Picture  of  the  Carl  Vinson  aircraft  carrier  [From  Ref  30] 


Figure  6:  The  aircraft  carrier  three-dimensional  model  in  four  view  angles. 


The  three-dimensional  wireframe  model  for  the  destroyer  was  based  on  a  1:1500 
scaled  drawing  of  the  Oscar  Austin  destroyer  (Arleigh  Burke  Class)  [Ref.  30].  This 
drawing  is  shown  in  Figure  7  and  a  picture  of  the  Oscar  Austin  destroyer  is  shown  in 
Figure  8.  The  destroyer  model  was  formed  of  92  vertices  and  57  planes  (see  the 
“destroyer.m”  program  in  Appendix  A).  Figure  9  shows  the  model  from  four  view  angles. 
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Figure  7:  Scaled  drawing  of  the  Oscar  Austin  destroyer  [From  Ref  30]. 
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Figure  8:  Picture  of  the  Oscar  Austin  destroyer  [From  Ref.30] 


Figure  9:  The  destroyer  three-dimensional  model  in  four  view  angles 


The  three-dimensional  wireframe  model  for  the  frigate  was  based  on  a  1:1200 
scaled  drawing  of  the  Rentz  frigate  (Oliver  Hazard  Perry  Class)  [Ref  30],  This  drawing 
is  shown  in  Figure  10  and  a  picture  of  the  Oscar  Austin  destroyer  is  shown  in  Figure  11. 
The  frigate  model  was  formed  of  130  vertices  and  66  planes  (see  the  “frigate.m”  program 
in  Appendix  A).  Figure  12  shows  the  model  from  four  different  view  angles. 


I 


,j£Ol„ 


‘O 


.  ^ . 

i'  4 


; . l^J. 


Figure  10:  Scaled  drawing  of  the  Rentz  frigate  [From  Ref.30]. 


Figure  1 1 :  Pictureiof  the  Rentz  frigate  [From  Ref  30]. 
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Figure  12:  The  frigate  three-dimensional  model  in  four  view  angles. 


The  three-dimensional  wireframe  model  for  the  merchant  Ship  was  based  on  a 
scaled  drawing  of  the  Sea  Isle  City  U.S.  tanker  [Ref.  31].  This  drawing  is  shown  in 
Figure  13  and  a  picture  of  the  Sea  Isle  City  U.S.  tanker  is  shown  in  Figure  14.  The 
merchant  model  was  formed  of  100  vertices  and  58  planes  (see  the  “merchant.m” 
program  in  Appendix  A).  Figure  15  shows  the  model  viewed  from  four  view  angles. 


Figure  13:  Scaled  drawing  of  the  Sea  Isle  City  tanker  [From  Ref.31]. 
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The  three-dimensional  wireframe  model  for  the  research  ship  was  based  on 
general  specifications  and  dimensions  extracted  from  RA^  Point  Sur  Cruise  Planning 
Manual  (see  Figure  16).  The  research  ship  model  was  formed  of  76  vertices  and  32 
planes  (see  the  “poitsur.m”  program  in  Appendix  A).  Figure  17  shows  the  model  viewed 
from  four  view  angles. 


Figure  16:  Picture  of  the  RA^  Point  Sur. 


Figure  17:  The  Point  Sur  three-dimensional  model  in  four  view  angles 


2.  Viewpoint  Control 


To  extract  silhouettes,  the  orientation  of  the  three-dimensional  ship  model  must 
be  specified.  It  was  possible  to  specify  the  viewpoint  with  the  MATLAB  “view” 
command  by  defining  azimuth  and  elevation  with  respect  to  the  axis  origin.  Azimuth  is  a 
polar  angle  in  the  x-y  plane,  with  positive  angles  indicating  counter-clockwise  rotation  of 
the  viewpoint.  Elevation  is  the  angle  above  (positive  angle)  or  below  (negative  angle)  the 
x-y  plane.  The  counter-clockwise  concept  for  the  azimuth  was  adopted  because  the 
viewing  azimuth  is  the  negative  of  the  ship’s  heading.  Therefore,  a  ship  with  heading  of 
30°  clockwise  is  equivalent  to  viewing  that  ship  with  azimuth  of  30°  counter-clockwise. 
The  diagram  in  Figure  1 8  illustrates  the  coordinate  system.  The  arrows  indicate  positive 
directions.  The  origin  was  assumed  to  be  located  approximately  in  the  center  of  gravity  of 
the  ship  model.  Only  the  portion  above  sea  level  was  considered.  Using  this  coordinate 
system,  we  can  verify  that  the  broadside  view  of  any  ship  model  corresponds  to  0°  in 
azimuth  and  0°  in  elevation,  in  this  situation  the  bow  direction  will  be  to  the  right 
(positive  x). 


Figure  18:  Diagram  illustrating  the  coordinate  system  and  the  ship  model  origin. 

One  view  of  each  of  the  five  modeled  ships  is  illustrated  in  Figure  19.  In  this 
figure,  the  azimuth  angle  is  -37.5  degrees  and  the  elevation  angle  is  30  degrees. 
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Figure  19:  Views  of  the  three-dimensional  ship  models  (az .=-37.5°,  elev.=30°). 


3.  Orthographic  Projection 

Once  the  three-dimensional  model  is  created  and  the  aspect  angle  is  set  using  the 
“view”  command,  a  silhouette  can  be  created  by  projecting  the  three-dimensional  ship 
model.  The  orthographic  method  projects  the  viewing  volume  as  a  rectangular 
parallelepiped  onto  a  plane,  i.e.,  relative  distance  from  the  camera  does  not  affect  the  size 
of  objects.  Using  orthographic  projection  it  is  possible  to  get  a  good  approximation  of  the 
real  process  of  image  generation  when  the  distance  from  the  object  to  the  camera  is  much 
greater  than  the  relative  deep  of  the  object  structural  points.  This  applies  to  our  task 
because  ships  are  generally  far  away  from  the  FLIR  sensors.  Figure  20  shows  some 
silhouettes  created  using  the  orthographic  projection  applied  to  the  images  shown  in 
Figure  19. 


Figure  20:  Silhouettes  created  using  the  images  of  Figure  19. 
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THE  REAL  FLIR  IMAGES  DATABASE 


1.  Domain  Issues 

We  also  obtained  real  images  taken  at  sea  using  the  AN/AAS-44V  Forward- 
Looking  Infrared  (FLIR)  sensor,  mounted  on  a  springboard  at  the  nose  of  the  SH-60B 
Rapid  Deployment  Kit  equipped  helicopter.  The  FLIR  images  were  available  through  a 
VHS-format  videotape  showing  several  ships.  However,  only  images  of  our  modeled 
ships  were  considered  for  our  analysis.  The  FLIR  images  show  good  contrast  and  were 
displayed  in  black-hot  (higher  temperatures  areas  are  black)  format.  Image  frames  were 
acquired  using  a  commercial  video-grabber  board  installed  in  a  PC-type  desktop 
computer. 

The  real  FLIR  images  that  we  extracted  from  the  videotapes  were  used  to  test  our 
recognition  system,  previously  trained  with  the  three-dimensional  ship  model  data.  This 
was  considered  particularly  important  since  presenting  the  classifier  with  new  images 
with  blurriness  and  an  unknown  target  viewpoint  is  challenging. 

Only  25  real  FLIR  images  were  used  for  testing  due  to  the  small  number  of 
modeled  ships  in  the  FLIR  tape.  These  were;  two  destroyer  images  (Figure  21),  four 
aircraft  carrier  images  (Figure  22),  15  merchant  ship  images  (Figure  23),  four  research 
ship  images  (Figure  24),  and  no  fiigate  images. 


Figure  21 :  Real  FLIR  images  of  destroyers 
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Figure  22:  Real  FLIR  images  of  aircraft  carrier  class. 


26 


Figure  24:  Real  FLIR  images  of  research  ship  class. 

The  quality  of  the  images  was  not  ideal  because  the  FLIR  system  projects 
alphanumeric  data  and  targeting  aids  onto  the  screen.  Figure  25  shows  one  of  the  FLIR 
images.  Specifically,  the  crosshairs  partially  obscure  the  ship  image  and  interfere  with  the 
classification  process.  A  segmentation  was  necessary  to  eliminate  the  backgroimd 
including  the  alphanumeric  data. 
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Figure  25:  FLIR  image  of  the  Arleigh  Burke  Destroyer. 


2.  Segmentation 

Prior  to  computing  the  moment  invariants  of  the  ship  of  the  Figure  25,  wq  must 
suppress  the  background  and  extract  the  ship  silhouette.  In  this  segmentation,  we 
employed  histogram  and  thresholding  techniques. 

We  assume  that  the  extracted  320x240  pixels  image  contains  one  ship  only.  This 
image  includes  the  ship,  water,  alphanumeric  data,  and  may  include  the  sky.  As  the  first 
step,  we  generated  the  gray-level  histogram  of  the  image  and  selected  a  threshold  level 
that  best  extracted  the  ship  from  the  water  region.  Figure  26  shows  the  histogram  of 
Figure  25. 
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In  this  image,  the  average  brightness  for  the  water  region  is  greater  than  that  for 
the  ship  region.  If  we  had  a  sky  region  in  this  picture,  the  values  in  the  sky  region  would 
be  lower  than  in  the  ship  region.  The  histogram  profile  was  analyzed  and  it  was  verified 
that  the  highest  peak  was  related  to  the  water  portion  of  the  image.  It  was  also  verified 
that  the  first  peak  left  of  the  highest  was  related  to  the  ship.  The  region  between  the  “ship 
peak”  and  the  “water  peak”  was  a  transition  region;  the  threshold  value  selected 
corresponded  to  the  minimum  in  the  transition  region.  In  Figure  26,  those  values  are 
water  peak=  95,  ship  peak=  65  and  threshold=  72.  The  original  image  was  thresholded  at 
that  value  and  a  binary  image  generated.  Figure  27  shows  the  result  for  Figure  25. 
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Figure  27:  Binary  image  thresholded  at  level  72  using  Figure  25. 

In  a  second  step,  we  eliminated  spurious  pixels  by  extracting  the  greatest 
connected  region  and  filling  the  holes  (see  segmentation.m  in  Appendix  A).  The  final 
ship  silhouette  found  for  Figure  25  is  shown  if  Figure  28. 


Figure  28:  Silhouette  found  after  cleaning  up  Figure  27. 
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D.  TRAINING  PHASE  OF  THE  NEURAL  NETWORK  CLASSIFIER 


The  neural  network  ship  classifier  required  a  training  phase  using  representative 
projective  views.  As  explained  in  Section  B  of  Chapter  n,  the  moment  functions 
invariance  under  reflection  helps  to  simplify  the  range  of  all  representative  views  distinct 
ship  views.  Three-dimensional  objects  which  possess  symmetry  about  a  plane,  such  as  a 
ship,  can  have  its  significant  range  of  distinct  views  for  azimuth  restricted  to  [-90°,  90°], 
where  0°  corresponds  to  the  broadside  ship  silhouette  [Ref  16].  Elevation  angles  were 
restricted  to  the  upper  hemisphere  with  45  degrees  as  the  upper  operational  limit,  as  the 
helicopter  will  be  viewing  the  ship  with  lower  elevation  angles. 

One  training  set  and  two  testing  sets  of  projections  were  generated.  The  training 
set  was  48  views  of  each  of  the  five  ship  types  taken  at  viewpoints  separated  by  15°  in  a 
180°x45°  sector;  i.e.,  { (0,  a),  0  =  -90°,  -75°,  -60°, ...,  75°;  a  =  0°,  15°,  30°,  45°  }  where  0 
and  a  represent  azimuth  and  elevation  angles  respectively.  Examples  of  the  training 
images  are  in  Figure  29  below.  In  this  figure,  the  elevation  angle  is  15°  and  the  azimuth 
angles  from  left  to  right  are:  -90°,  -60°,  -30°,  0°,  30°,  and  60°.  The  ship  types  from  top  to 
bottom  are  aircraft  carrier,  destroyer,  fngate,  merchant,  and  research  ship. 
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Figure  29:  Examples  of  the  training  images  for  each  modeled  ship 


E.  TESTING  PHASE 

There  were  two  testing  sets;  41400  silhouettes  projected  from  the  three- 
dimensional  ship  models  and  the  25  real  FLIR  images.  The  first  test  set  contains  8280 
views  of  each  of  the  five  ship  models  taken  at  viewpoints  separated  by  1°  in  azimuth  and 
r  in  elevation;  i.e.,  {(e,a),  0  =  -90°,  -89°,  -88°, ...,  89° ;  a  =  0°,  1°,  2°,  ...,  45° }  where  0 
and  a  represent  azimuth  and  elevation  angles  respectively.  Although  this  set  contains  the 
training  set,  the  number  of  training  views  was  very  small  (96  views)  to  compromise  the 
simulation  results.  The  real  image  test  set  was  described  earlier. 


F. 


PROGRAMS  DEVELOPED 


All  16  programs  were  written  in  MATLAB.  They  can  be  divided  in  four 
categories  (see  Table  1  and  Appendix); 

-  Three-dimensional  ship  modeling:  specialized  functions  used  to  create 
the  three-dimensional  models  used  in  other  programs; 

-  Moment  invariant  computation:  programs  used  to  calculate  the 
moment  invariants  of  a  specific  ship  silhouette; 

-  Neural  network  training:  programs  used  to  train  the  neural  network 
implemented;  and 

Testing:  programs  used  to  evaluate  the  performance  of  the  system 
implemented. 
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Name 

Type 

Function 

aircarrier.m 

Three- 

Create  the  model  for  aircraft  carrier 

destroyer.m 

dimensional 

ship 

Create  the  model  for  destroyer 

frigate.m 

modeling 

Create  the  model  for  frigate 

pointsur.m 

programs 

Create  the  model  for  research  ship 

merchant.m 

Create  the  model  for  merchant  ship 

findInputMomSet.m 

Moment 

Returns  the  12-element  input  set 

find_mom_functions.m 

invariants 

computation 

Returns  the  six  moment  functions  values 

find_moment.m 

programs 

Returns  the  central  moment 

find_centroid.m 

Returns  the  centroid  of  a  silhouette 

mainShipRecon.m 

Neural 

network 

Creates  and  trains  a  neural  network 
responsible  for  recognizing  ship  types 

findInputSet.m 

training 

programs 

Returns  all  the  silhouettes  to  be  used  by 
the  neural  network  during  training  phase 

interface.m 

Testing 

programs 

Create  a  Graphical  User  Interface  to 
evaluate  the  system  implemented.  Three 
ship  silhouettes  are  shown  in  the 
interface:  (l)the  original  silhouette,  (2)the 
rotated,  scaled  and  noisy  silhouette 
defined  by  the  user,  and  (3)the  neural 
network  guessed  silhouette 

plotSilhouette.m 

Draws  the  ship  silhouette  inside  the 
Graphical  User  Interface 

segmentation.m 

Segments  a  FLIR  real  image  using  a 
histogram  and  threshold  technique 

createTestSet.m 

Creates  05  mat  files  containing  the 
silhouettes  of  each  ship  for  increments  of 
one  degree  in  azimuth  and  elevation, 
then  plots  the  errors 

findResultSet.m 

Returns  a  vector  with  the  size  of  all  the 
viewangles  being  tested,  where  ”1"  will 
mean  misclassified  and  "0"  will  mean 
correct  classified 

Table  1:  Programs  implemented 
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V.  RESULTS  FROM  EXPERIMENTATION 


A.  EVALUATION  STRUCTURE 

Experiments  were  carried  out  in  order  to  evaluate  the  proposed  system.  For  this 
purpose,  the  five  ship  models  detailed  in  chapter  IV  were  used.  The  parameters  generally 
used  to  characterize  the  overall  performance  of  an  automatic  target  recognition  system 
are  the  probability  a  given  ship  is  correctly  recognized  (recall)  and  probability  a  given 
type  identification  is  correct  (precision).  By  classification  is  meant  recognition,  the 
determination  of  the  target  type.  Performance  data  is  generally  given  in  the  form  of  a 
confusion  matrix  together  with  the  size  of  the  feature  vector  and  database. 

Using  the  probability  of  correct  classification  as  a  reference  and  the  simulated 
three-dimensional  ship  model  images  as  a  database,  it  was  possible  to  iteratively  optimize 
our  system.  This  optimization  was  achieved  in  three  distinct  experiments. 

B.  FIRST  EXPERIMENT 

The  first  experiment  was  implemented  as  described  in  chapter  IV.  Experiments 
were  performed  with  simulated  12-components  moment  invariants  signatures  from  the 
models.  The  six  moment  invariants  of  the  solid  silhouette  and  the  six  moment  invariants 
of  the  boundary  make  up  the  signature  vector,  as  described  in  chapter  n.  The  training 
signatures  were  generated  from  images  of  the  three-dimensional  models  taken  at  a  regular 
pace  of  15  degrees  in  azimuth  increments  and  in  four  different  viewing  elevation  angles: 
0°,  15®,  30°  and  45°.  The  total  training  set  was  48  images. 

The  testing  set  contained  the  images  of  the  three-dimensional  models  taken  at  a 
regular  pace  of  1°  in  azimuth  increments  and  using  the  same  four  elevation  angles  used 
during  training.  The  neural  network  used  was  small,  with  only  20  hidden  neurons  and  a 
single  hidden  layer.  Backpropagation  was  used  as  the  training  technique.  The  network 
was  successfully  trained.  This  network  yielded  90.1%  discrimination  leaving  an  overall 
approximate  10%  error  rate.  Details  of  the  experiment  are  summarized  in  Table  2. 


37 


Network  Parameters 

Network- 1 

#  input  nodes 

12 

#  nodes  in  hidden  layer 

20 

#  output  nodes 

5 

#  training  set 

240  (12  X  4  X  5) 

#  test  set 

3600(180x4x5) 

Accuracy 

90.1% 

Table  2:  Neural  network  for  experiment  1 


The  first  experiment  enabled  us  to  analyze  the  azimuth  behavior  of  the  neural  net 
generalization  capability.  For  this  purpose,  we  plotted  for  each  ship  type  the  classification 
error  percentage  with  respect  to  azimuth  (Figures  30  to  32). 


Aircrat  Carrier  Testing  Results  (adding  all  4  elevation  angles) 


azimuth  angle 


Figure  30:  Accuracy  with  respect  to  azimuth  for  aircraft  carrier 
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Pointsur  Testing  Results  (adding  all  4  elevation  angles) 
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Merchant  Testing  Results  (adding  all  4  elevation  angles) 
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Figure  32:  Accuracy  with  respect  to  azimuth  variation  for  Point  Sur  and  merchant 


Analyzing  the  results  of  these  graphs,  we  can  verify  that  the  major  errors  were 
related  to  high  azimuth  angles.  This  error  pattern  helped  us  to  decide  to  use  an  unevenly 
spaced  training  set.  In  the  second  experiment,  we  increased  the  number  of  azimuth 
training  angles  by  using  small  steps  (5°)  in  high  azimuth  angles. 

C.  SECOND  EXPERIMENT 

In  the  second  experiment,  we  increased  the  number  of  training  azimuth  angles  by 
using  an  unevenly  spaced  training  set.  For  low  azimuth  angles  (-45°  to  45°),  we  kept  the 
15°  step,  but  for  high  azimuth  angles  we  chose  a  smaller  step  of  5°  .The  new  training  set 
contained  96  views  of  each  of  the  five  ship  types  taken  at  the  following  viewpoints:  {(0, 
a),  0  =  -90,  -85,  -80,  -75,  -70,  -65,  -60,  -55,  -50,  -45,  -30,  -15,  0,  15,  30,  45,  50,  55,  60, 
65,  70,  75,  80,  85;  a  =  0,  15,  30,  45},  where  0  and  a  represent  azimuth  and  elevation 
angles  respectively. 

The  neural  network  and  testing  set  were  the  same  used  in  the  first  experiment.  We 
obtained  91.2%  accuracy.  We  also  tried  another  neural  network  architecture  with  30 
hidden  nodes  and  the  classification  rate  improved  to  94.8%.  Details  of  these  experiments 
are  summarized  in  Table  3. 


Network  Parameters 

Network-2 

Network-3 

#  input  nodes 

12 

12 

#  nodes  in  hidden  layer 

20 

30 

#  output  nodes 

5 

5 

#  training  set 

480  (24  X  4  X  5) 

480  (24  X  4  X  5) 

#  test  set 

3600(180x4x5) 

3600(180x4x5) 

Accuracy 

91.2% 

94.8% 

Table  3:  Neural  networks  for  experiment  2 


As  we  can  see  in  Table  3,  with  a  bigger  training  set  and  30  hidden  neurons  the 
neural  network-3  )delded  the  best  accuracy.  Figures  33  to  35  show  the  classification 
error  percentage  with  respect  to  azimuth  for  network-3. 
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Aircrat  Carrier  Testing  Results  (30  hidden  nodes  and  adding  all  4  elevation  angles) 
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Destroyer  Testing  Results  (30  hidden  nodes  and  adding  all  4  eIe\a(tion  angles) 
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Figure  33:  Accuracy  with  respect  to  azimuth  for  aircraft  carrier  and  destroyer 
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Frigafte  Testing  Results  (30  hidden  nodes  and  adding  all  4  elesation  angles) 
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Pointsur  Testing  Results  (30  hidden  nodes  and  adding  all  4  ele\ation  angles) 
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Figure  34:  Accuracy  with  respect  to  azimuth  for  frigate  and  Point  Sur 


43 


Merchant  Testing  Results  (30  hidden  nodes  and  adding  all  4  elevation  angles) 


Figure  35:  Accuracy  with  respect  to  azimuth  for  merchant  ship 


To  analyze  the  elevation  behavior  of  the  neural  network-3,  we  graphed,  for  each 
ship  type,  the  classification  error  percentage  with  respect  to  elevation.  Views  of  the  three- 
dimensional  model  were  calculated  in  increments  of  1°  in  elevation,  starting  with  0°  and 
ending  with  45°,  and  using  the  same  24  azimuth  angles  used  during  training  phase  of 
network-3  (Figures  36  to  38). 

Analyzing  Figures  36  to  38,  we  can  verify  that  the  major  errors  were  related  to 
small  elevation  angles.  This  elevation  error  pattern  inspired  us  to  perform  the  third 
experiment,  this  time  we  increased  the  number  of  elevation  training  angles.  The  idea  was 
to  improve  the  classification  accuracy  of  network-3  by  training  with  more  elevation 
angles. 
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Frigate  Testing  Results(30  hidden  nodes  and  adding  all  24  azimuth  angles) 


P  ointsur  Testing  Results(30  hidden  nodes  and  adding  all  24  azim  uth  angles) 


Figure  37: Accuracy  with  respect  to  elevation  for  frigate  and  Point  Sur 
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Merchant  Testing  Results(30  hidden  nodes  and  adding  all  24  azimuth  angles) 


D.  THIRD  EXPERIMENT 

To  improve  the  classification  accuracy  of  network-3  we  increased  the  number  of 
elevation  angles  used  for  training  in  a  third  experiment.  The  new  training  set  contains  168 
views  of  each  of  the  five  ship  types  taken  at  the  following  viewpoints:  {(0,  a),  0  ==  -90°,  - 
85°,  -80°,  -75°,  -70°,  -65°,  -60°,  -55°,  -50°,  -45°,  -30°,  -15°,  0°,  15°,  30°,  45°,  50°,  55°,  60°, 
65°,  70°,  75°,  80°,  85°;  a  =  0°,  7°,  15°,  22°,  30°,  45°},  where  0  and  a  represent  azimuth 
and  elevation  angles  respectively.  The  neural  network  and  testing  set  were  the  same  used 
before.  We  obtained  85.4%  discrimination  leaving  an  overall  14.6%  error  rate  (Table  4). 

This  new  architecture  presented  worse  results  because  the  neural  network  did  not 
converge  during  the  training  phase  using  the  backpropagation  algorithm.  This 
convergence  problem  was  due  to  the  large  number  of  the  training  samples. 
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Network  Parameters 

Network-4 

#  input  nodes 

12 

#  nodess  in  hidden  layer 

30 

#  output  nodes 

5 

#  training  set 

720  (24  X  6  X  5) 

#  test  set 

3600  (180x4x5) 

Accuracy 

85.4% 

Table  4:  Neural  network  for  experiment  3 


E.  FOURTH  EXPERIMENT 

The  third  experiment  has  proved  that  the  neural  network-3  (accuracy  of  94.8%) 
was  our  optimized  solution.  Thus,  it  was  adopted  for  testing  with  the  whole  set  of  model 
images  and  with  the  real  FLIR  images. 

The  fourth  experiment  was  to  apply  our  approach  to  a  large  final  testing  set.  This 
final  set  contained  the  images  of  the  three-dimensional  ship  models  taken  at  a  regular 
pace  of  1°  in  azimuth  and  elevation  increments,  totaling  41400  images  (180  x  46  x  5). 
This  experiment  was  performed  using  the  neural  network-3,  from  the  second  experiment. 
The  average  classification  accuracy  achieved  was  87.3%.  Table  5  shows  the  confusion 
matrix  of  counts. 
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X  Inferred 
?\Type 
InpuK 
Type.  \ 

Aircraft 

Carrier 

Destroyer 

Frigate 

Point 

Sur 

Merchant 

Accuracy 

Aircraft 

Carrier 

6711 

201 

425 

393 

550 

81.1% 

Destroyer 

318 

7301 

397 

257 

7 

88.2% 

Frigate 

345 

788 

6809 

217 

121 

82.2% 

Pointsur 

67 

146 

177 

7873 

17 

95.1% 

Merchant 

297 

188 

291 

49 

7455 

90.0% 

Precision 

86.7% 

84.7% 

84.1% 

89.6% 

91.5% 

Overall  probability  of  classification: 

87.3% 

Table  5:  Confusion  matrix  for  41,400  views  of  modeled  images  and  network-3 


F.  FIFTH  EXPERIMENT 

We  also  ran  the  system  on  the  25  real  FLIR  images,  as  described  in  chapter  IV. 
This  experiment  was  also  performed  using  the  neural  network-3,  from  the  second 
experiment.  The  average  classification  accuracy  achieved  was  68%.  Table  6  shows  the 
confusion  matrix  of  counts. 


N^^Inferred 
.  \Type 
InpuNv 
Type  X 

Aircraft 

Carrier 

Destroyer 

Point 

Sur 

Merchant 

Recall 

Aircraft 

Carrier 

3 

0 

0 

1 

75% 

Destroyer 

0 

2 

0 

100% 

Pointsur 

0 

1 

3 

75% 

Merchant 

6 

0 

9 

60% 

Precision 

33.3% 

66.6% 

100% 

90% 

Overall  probability  of  classification: 

68% 

Table  6:  Confusion  matrix  for  real  ship  FLIR  images  and  network-3. 
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Figures  39  to  44  show  all  the  FLIR  images  and  respective  results.  The  first 
columns  show  the  original  FLIR  images;  the  middle  columns  show  the  silhouettes  after 
the  segmentation  process;  and  the  third  columns  show  the  neural  network’s  guess. 


Figure  39:  Classification  results  for  aircraft  carrier  FLIR  images 
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Figure  40:  Classification  results  for  merchant  FLIR  images 
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Figure  41 :  Classification  results  for  merchant  FLIR  images 
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Figure  42:  Classification  results  for  merchant  FLIR  images 
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Figure  43:  Classification  results  for  point  sur  FLIR  images 
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Figure  44;  Classification  results  for  destroyer  FLIR 
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VI.  CONCLUSIONS 


This  thesis  explored  a  moment-based  method  for  ship-type  recognition. 
Numerical  simulations  were  carried  out  with  a  set  of  five  three-dimensional  ship  models. 
Moment-invariant  signatures  were  used  as  the  input  to  a  neural-network  classifier.  The 
classifier  achieved  a  87.5%  correct  classification  rate  (within  the  set  of  test-models)  for  a 
complete  range  of  point  of  view  around  the  input  ship  model.  The  success  is  due  to  a 
combination  of  a  robust  feature  extraction  and  the  neural-network  generalization 
capability. 

A  test  with  25  real  FLER  ship  images  was  also  done.  Experimental  results  were 
worse  due  to  the  noisy  extracted  silhouettes.  The  maximum  classification  accuracy  of 
68%  should  only  be  considered  a  rough  approximation  to  the  sort  of  accuracies  one  can 
expect  from  a  fully  operational  classifier. 

A  larger  database  of  real  FLIR  ship  images  needs  to  be  tested.  The  acquisition 
process  should  eliminate  undesired  alphanumeric  data  superimposed  on  the  FLIR  images. 
Using  this  database,  our  system  could  be  trained  using  the  real  FLIR  images  from 
different  viewpoints,  and  this  could  provide  better  performance  on  new  real  images. 
Another  recommendation  is  to  investigate  different  segmentation  algorithms  capable  of 
addressing  some  predictable  distortions  like  the  ship  shadow/reflection  on  the  sea  surface 
and  the  smoke  coming  out  from  stacks.  As  we  can  see  in  Figure  35,  the  merchant  FLIR 
images  presented  the  worst  accuracy  due  to  these  distortions. 

In  conclusion,  this  thesis  demonstrated  how  ship  recognition  using  models  is 
complicated  by  the  imaging  process,  which  involves  a  viewpoint-dependent  two- 
dimensional  projection  of  three-dimensional  ship  model.  As  a  consequence,  the 
appearance  of  a  ship  in  an  image  can  vary  greatly  with  its  aspect  and  scale.  Ship 
recognition  from  infrared  images  is  further  complicated  because  the  extracted  ship 
silhouettes  can  be  noisy  due  to  the  distortions  caused  by  shadows,  smoke,  and  other 
factors. 


This  thesis  demonstrated  the  potential  of  a  simple  algorithm  for  this  particular 
application.  The  modest  requirements  in  terms  of  computer  and  FLIR  hardware  of  this 
system  show  great  potential  for  providing  a  recognition  system  to  a  variety  of  users. 
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APPENDIX  A.  PROGRAM  LISTING 


Naval  Postcraduate  Schccl 


^  Type  :  Function 

'<?  Narne  :  aircarmen.i'n 

%  Function  ;  returns  the  vertices  and  faces  to  be  used  by  the  M^TLAB  functio 
"patch"  and  construct  the  aircraft  carrier  3-D  model 
%  Date  G1  march  2  0 Cl 

%  Version  :  1.0 

%  Author  :  Jorge  Amaral  Alves,  LCDR  (Brazilian  Navy) 


function  [verts , faces]  =aircarrier ( )  ; 


verts= [122 , 

-7, 

0  %1 

122, 

7, 

0  %2 

116, 

1, 

-7  %3 

114, 

0, 

-12 

116, 

-1, 

-7  %5 

77, 

-10, 

0  %6 

77, 

-10, 

-7  %7 

77, 

-8, 

-12  %S 

68, 

-10, 

-7  %9 

68, 

-21, 

0  %1G 

-23, 

-21, 

0  %11 

-23, 

-10, 

-7  %i2 

-62, 

-10, 

0  %13 

-62, 

-10, 

-7  %14 

-60, 

-8, 

-12  % 

-62, 

10, 

0  %i6 

-62, 

10, 

-7  %17 

-60, 

8, 

-12  % 

-36, 

12, 

0 

-36, 

10, 

-7  %20 

-28, 

21, 

0  %21 

-28, 

10, 

-7  322 

38, 

21, 

0  %23 

38, 

10, 

-7 

65, 

24, 

0 

65, 

10, 

-  7  %2 6 

77, 

10, 

0  ^27 

77, 

10, 

-7  %23 

77, 

8, 

-12  %29 

-7, 

-19, 

0  %30 

-7, 

-19, 

9  %3i 

-7, 

-15, 

9  %32 

-7, 

-15, 

0  %3  3 

7, 

-19, 

0  %34 

7, 

-19, 

9  %35 

7, 

-12, 

9  %36 

7, 

-12, 

0  %37 

-1, 

-17, 

9  %3S 

-1, 

-17, 

20  %39 

-1, 

-15, 

20  %40 

-1, 

-15, 

9  %41 

1/ 

-17, 

9  %42 

59 


1, 

-17, 

20  %4 

3 

1, 

-15, 

9  %4 

-1 

1, 

-15, 

20]  ;% 

4  5 

%nuLnber 

- 

32  12 

planel  = 

[1, 

2,3,4, 

5,5 

,5,5,5 

,5,5,5] ; 

plane2  = 

[1, 

5,S,S, 

S,6 

,5,6,6 

,6,6,6] ; 

planes  = 

[5, 

^,1 .1 , 

1.1 

,  7, ?, ? 

,7,?,?]; 

plane4  = 

[4, 

5,?,S, 

S,S 

,8,8,8 

,8,8,8] ; 

planes  = 

[6, 

9,  ?, ?, 

1.1 

,7,?,? 

,7,?,?] ; 

planes  = 

[6, 

10,S,S 

,9, 

9,9,3, 

9, 3, 9, 9]  ; 

plane?  = 

[9, 

10,11, 

12, 

L2,12, 

12,12,12, 

12,12,: 

12]  ; 

planes  = 

[11 

,12,13 

,13 

,13,13 

,13,13,13 

,13 

,13 

,13] 

; 

planes  = 

[12 

/13,14 

,14 

,14,14 

r 14, 14, 14 

,14 

14 

,14] 

/ 

planelO 

= 

[7,9,12,1 

4,15,8 

,8, 8, 8, 8, 

8,8 

]  ; 

planell 

= 

[13,16 

,17 

,14,14 

,14,14,14 

,14 

,14 

,  14, 

14] 

planel2 

[14,17 

,  18 

,15,15 

,15,15,15 

,15 

,15 

,15, 

15] 

planelS 

= 

[16,19 

,20 

,17,17 

,17,17,17 

,17 

,17 

,17, 

17] 

planel4 

= 

[19,21 

,22 

,20,20 

,20,20,20 

,20 

,20 

,20, 

20] 

planelS 

= 

[21,23 

,24 

,22,22 

,22,22,22 

,22 

,22 

,22, 

22] 

planelS 

= 

[23,25 

,2S 

,24,24 

,24,24,24 

,24 

,24 

,24, 

24] 

planel? 

= 

[25,27 

,2S 

,26,26 

,26,26,26 

,26 

,26 

,26, 

26] 

planelS 

[26,27 

,28 

,28,28 

,28,28,28 

,28 

,28 

,28, 

28] 

planelS 

[27,3, 

28, 

28,28,; 

28,28,28,: 

28, 

28, 

28,28]  ; 

plane20 

= 

[27,2, 

3,3 

,3,3,3 

,3, 3, 3, 3,: 

3]  ; 

plane21 

[28,3, 

4,2 

9,23,23,23,23,2 

9,2 

9,29,23 

]  ; 

plane22 

= 

[17,20 

,22 

,24,26 

,28,29,18 

,  18 

,18 

,18, 

18] 

f lightdeck 

=  [1, 

6,10,11,13,16,13,21,23,2 

5,27,2] 

bottom  = 

[4, 

8,15,18,2S 

),23,2S 

), 29, 29,29,25 

1,29]  ; 

towerl  = 

[30 

,31,32 

33, 

33,33, 

33,33,33, 

33, 

33, 

33] 

; 

tower 2  = 

[30 

,31,35 

34, 

34,34, 

34,34,34, 

34, 

34, 

34] 

/ 

towers 

= 

[34,35 

,36 

,37,37 

,37,37,37 

,37 

,37 

,37, 

37] 

tower 4  = 

[32 

,33,3? 

36, 

36,36, 

36,36,36, 

36, 

36, 

36] 

t 

towers  = 

[31 

,32,36 

35, 

35,35, 

35,35,35, 

35, 

35, 

35] 

, 

mastl  = 

[3S 

,33,40 

41, 

41,41, 

41,41,41, 

41, 

41, 

41] 

; 

mast2  = 

[3S 

,33,43, 

42, 

42,42, 

42,42,42, 

42, 

42, 

42] 

, 

masts  = 

[42 

,44,45, 

43, 

43,43, 

43,43,43, 

43  , 

43, 

43] 

, 

mast4  = 

[41 

,44,45, 

40, 

40,40, 

40,40,40, 

40, 

40, 

40] 

; 

masts  = 

[33 

,43,45, 

40, 

40,40, 

40,40,40, 

40, 

40, 

40] 

; 

%faces 

faces= [  planel 

plane2 

planes 

plane4 

planes 

plane6 

plane? 

planes 

planes 

planelO 

planell 

planel2 

planelS 

planel4 

planelS 

planelG 

planel? 

planelS 
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planel9 

plane20 

plane21 

plane22 

f lightdeck 

bottom 

towerl 

tower2 

towers 

tower4 

towers 

mastl 

mast2 

masts 

mast4 

masts  ]  ; 


End  of  file  aircarrier.m 


61 


%  Nav.?*i  Postgraduate  School  -  CA 

:  Function 
:  dsstroyer.m 

>n  :  returns  the  vertices  and  faces  to  be  used  b\ 
’’patch"  in  order  to  construct  the  destroyer  3 
%  Date  01  march  2001 
%  Version  :  1.0 

%  Author  :  Jorge  Amaral  Alves,  LCBR  (Brazilian  Navy) 


function  [verts, faces] =destroyer { ) ; 


=  [140 

,  0, 

5  %1 

90, 

-16 

,  3  %2 

90, 

-16 

,  -10  %3 

125 

,  0, 

-10  %4 

50, 

-20, 

0  %5 

50, 

-20, 

-10  %6 

44, 

-20, 

0  %7 

44, 

-20, 

5  %8 

00 

ro 

-20, 

5  %9 

29, 

-20, 

5  %10 

8,  - 

20, 

5  %11 

8,  - 

20,  0%12 

-8, 

-20, 

0  %13 

“8, 

-20, 

-10  %14 

‘8, 

-20, 

5  %15 

-55, 

-14, 

5  %16 

-55, 

-14, 

-3  %17 

-55, 

-13 

-10  %1 

-85, 

-10, 

-3  %19 

-84, 

-9, 

-10  %20 

-85, 

10, 

-3  %21 

-84, 

9, 

-10  %22 

-55, 

14, 

-3  %23 

-55, 

13, 

-10  %24 

-55, 

14, 

5  %25 

-8, 

20, 

5  %26 

-8, 

20, 

0  %27 

-8, 

20, 

-10  %28 

8, 

20, 

0  %2S 

8, 

20, 

5  %30 

29, 

20, 

5  %31 

38, 

20, 

5  %32 

44, 

20, 

5  %33 

44, 

20, 

0  %34 

50, 

20, 

0  %35 

50, 

20, 

-10  %36 

90, 

16, 

3  %37 

90, 

15, 

-10  %3S 

50, 

-11, 

0  %39 

50, 

11, 

0  %40 

50, 

-6, 

5  %41 

44, 

-6, 

5  %42 

44, 

-11, 

5  %43 

38, 

-20, 

5  %44 

44, 

-20, 

5%4  5 

50, 

-11, 

5%46 

T  the  MAT LAB  functio 
“D  model 


62 


50, 

6, 

5%4  7 

44, 

6, 

5746 

44  , 

11, 

5%49 

38, 

20, 

5%50 

44, 

20, 

5751 

50, 

11, 

5%52 

50, 

-6, 

14%53 

44  , 

-6, 

14%54 

44, 

6, 

14%55 

50, 

5, 

14%55 

38, 

-19, 

20%57 

44, 

-11, 

20%58 

44, 

11, 

20159 

38, 

19, 

20%60 

29, 

19, 

20%61 

25, 

11, 

20%62 

25, 

-11, 

20%63 

29, 

-19, 

20164 

0, 

-10, 

5%65 

20, 

-10, 

5%66 

20, 

10, 

5167 

0, 

10, 

5  %  6  8 

24, 

11, 

5169 

24, 

-11, 

5%70 

6, 

-4, 

20%71 

14, 

-4, 

20172 

14  , 

4, 

20%73 

6,  4,  20%74 


0, 

-10, 

0175 

0, 

10, 

0%76 

-32, 

-10, 

5177 

-8, 

-10, 

517S 

-8, 

10, 

5%79 

-32, 

10, 

57  80 

-22, 

-4, 

20%81 

-14, 

-4, 

20182 

-14  , 

4, 

20183 

-22, 

4, 

20%84 

29, 

-2, 

20135 

33, 

-2, 

20%86 

33, 

2, 

20187 

29, 

2, 

20186 

25, 

-2, 

50189 

29, 

-2, 

50190 

29, 

2, 

50191 

25, 

2, 

50] ;%9 

%riumber  =  3  2  12 

planel  =  [1, 2, 3, 4, 4, 4, 4, 4, 4,4,4]  ; 

plane2  =  [2, 5, 6, 3, 3, 3, 3, 3, 3,3,3]  ; 

plane3  =  [5,7,8,9,10,11,12,13,14,6,6] ; 

plane4  =  [13,15,16,17,18,14,14,14,14,14,14] 

plane5  =  [17,19,20,18,18,18,18,18,18,18,18] 

plane6  =  [19,21,22,20,20,20,20,20,20,20,20] 

plane7  =  [21,23,24,22,22,22,22,22,22,22,22] 

plane8  =  [23,25,26,27,28,24,24,24,24,24,24] 

plane9  =  [27,29,30,31,32,33,34,35,36,28,28] 
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planelO  =  [35,37,38,36,36,36,36,36,36,36,36] 

planell  =  [37,1,4,38,38,38,38,38,38,38,38]; 

botton  =  [4,3,6,14,18,20,22,24,28,36,38] ; 

deckl  =  [1,2,37,37,37,37,37,37,37,37,37]; 

deck2  =  [2,5,35,37,37,37,37,37,37,37,37]; 

deck3  =  [5,7,39,39,39,39,39,39,39,39,39]; 

deck4  =  [34,35,40,40,40,40,40,40,40,40,40]; 

decks  =  [46,41,42,43,9,8,8,8,8,8,8]; 

deck6  =  [47,48,49,32,51,52,52,52,52,52,52]; 

deck7  =  [53,54,55,56,56,56,56,56,56,56,56]; 

decks  =  [57,58,59,60,61,62,63,64,64,64,64]; 

deck9  =  [10,11,65,66,67,68,30,31,69,70,70]; 

decklO  =  [13,12,75,76,29,27,27,27,27,27,27] 

deckll  =  [16,15,78,77,80,79,26,25,25,25,25] 

deckl2  =  [19,17,23,21,21,21,21,21,21,21,21] 

laterall  =  [8,46,39,7,7,7,7,7,7,7,7],- 

lateral2  =  [39,46,41,53,56,47,52,40,40,40,40] 

laterals  =  [40,52,33,34,34,34,34,34,34,34,34] 

lateral4  =  [54,53,41,42,42,42,42,42,42,42,42] 

laterals  =  [55,56,47,48,48,48,48,48,48,48,48] 

laterals  =  [58,59,49,48,55,54,42,43,43,43,43] 

lateral7  =  [57,58,43,9,9,9,9,9,9,9,9]; 

laterals  =  [60,59,49,32,32,32,32,32,32,32,32] 

laterals  =  [64,57,9,10,10,10,10,10,10,10,10],- 

laterall0=  [61,60,32,31,31,31,31,31,31,31,31] ; 
lateralll=  [63, 64, 10, 70, 70, 70, 70, 70, 70, 70, 70]  ,- 
laterall2=  [62,61,31,69,69,69,69,69,69,69,69] 
laterall3=  [62,63,70,69,69,69,69,69,69,69,69] 
laterall4=  [65,11,12,75,75,75,75,75,75,75,75] 
laterall5=  [68,30,29,76,76,76,76,76,76,76,76] 
laterall6=  [65,75,76,68,68,68,68,68,68,68,68] 
laterall7=  [15,13,27,26,26,26,26,26,26,26,26] 
laterall8=  [16,17,23,25,25,25,25,25,25,25,25]  ; 
stackl  =  [71,72,66,65,65,65,65,65,65,65,65],- 

Stack2  =  [74,73,67,68,68,68,68,68,68,68,68],- 

stack3  =  [65,71,74,68,68,68,68,68,68,68,68]; 

Stack4  =  [66,72,73,67,67,67,67,67,67,67,67]; 

Stacks  =  [71,72,73,74,74,74,74,74,74,74,74],- 

Stack6  =  [77,81,82,78,78,78,78,78,78,78,78]; 

Stack7  =  [80,84,83,79,79,79,79,79,79,79,79],- 

stackS  =  [77,81,84,80,80,80,80,80,80,80,80]; 

Stack9  =  [78,82,83,79,79,79,79,79,79,79,79],- 

stackio  =  [81,82,83,84,84,84,84,84,84,84,84]; 

mastl  =  [85,86,90,89,89,89,89,89,89,89,89],- 

mast2  =  [86,87,91,90,90,90,90,90,90,90,90],- 

mast3  =  [87,88,92,91,91,91,91,91,91,91,91],- 

mast4  =  [88,85,89,92,92,92,92,92,92,92,92],- 

mastS  =  [89,90,91,92,92,92,92,92,92,92,92],- 

%f aces 

f aces= [planel 
plane2 
planes 
plane4 
planes 
planes 
planeV 
planes 
plane9 


planelO 

planell 

botton 

deckl 

deck2 

decks 

deck4 

decks 

decks 

deck? 

decks 

decks 

decklO 

deckll 

deckl2 

laterall 

lateral2 

laterals 

lateral4 

laterals 

laterals 

lateral? 

laterals 

laterals 

laterallO 

lateralll 

lateralis 

lateralis 

laterall4 

lateralis 

lateralis 

laterall? 

lateralis 

stackl 

Stacks 

Stacks 

stack4 

Stacks 

Stacks 

Stack? 

Stacks 

Stacks 

stacklO 

mastl 

masts 

masts 

mast4 

masts]  ; 


%  End  of  file  destroys 


r ,  m 
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(A'J  oV 


Naval  Postgraduate  School  -  CA 

:  Function 
:  frigate. in 

on  :  returns  the  vertices  and  faces  to  be  used  by  the  M 
%  "patch"  in  order  to  construct  the  frigate  3-D  niodel 

%  Date  Q1  march  2001 
%  Version  :  l.C 

%  Author  :  Jorge  Amaral  Alves,  LCDR  (Brazilian  Navy) 


iction 


function  [verts , faces] ^frigate ( ) ; 


=  [128, 

0, 

16 

%1 

96, 

-13 

,  14 

%2 

96, 

-7, 

0 

%3 

112, 

0, 

0 

%4 

92, 

-15, 

12 

%5 

92, 

-9, 

0 

%6 

60, 

-20, 

10 

%7 

60, 

-18, 

,  0 

%8 

58, 

-20, 

10 

%9 

48, 

-20, 

10 

%10 

-26, 

-20, 

10 

%  1 1 

-60, 

-20, 

10 

%12 

-60, 

-18, 

0 

%13 

-102, 

-15, 

10 

%  1 4 

-98, 

-13, 

0 

%15 

-102, 

15, 

10 

%16 

-98, 

13, 

0 

%17 

-60, 

20, 

10 

%18 

-60, 

18, 

0 

%19 

-26, 

20, 

10 

%2G 

48, 

20, 

10 

%21 

58, 

20, 

10 

%22 

60, 

20, 

10 

%2; 

60,18 

,0 

%24 

92,15 

,12 

%2  5 

92,9, 

0  %26 

96,13 

,14 

%27 

96,7, 

0  %28 

60, -13,10 

%2  9 

60,13 

,10 

%3  0 

-12,- 

12,10 

%31 

8,-12 

,10 

%32 

8,-16 

,10 

%3  3 

-12,12,10 

%34 

8,12, 

10 

%35 

8,16, 

10 

%3  6 

58,-20,20 

%37 

48,-20,20 

%3S 

48,-16,20 

%39 

8,-16 

,20 

%40 

8,-12 

,20 

%41 

-12,- 

12,20 

%42 

-26,  - 

20,20 

%43 

-60,  - 

20,20 

%44 

-60,20,20  %45 

-26,20,20  %46 
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-12,12,20 

■in/ 

8,12,20 

%4  S 

8,16,20 

%49 

48,16,20 

i50 

48,20,20 

%  5 1 

58,20,20 

%52 

60,13,20 

%53 

60, -13,20 

%54 

57,-13,25 

%5  5 

43, -13,25 

%56 

43,13,25 

%57 

57,13,25 

%56 

60,7,25 

%5  9 

60,-7,25 

%6  0 

49,-20,20 

%6i 

49, -20,12 

%62 

51,-20,12 

%6  3 

51, -20,20 

%64 

52,-20,20 

%65 

52, -20,12 

%56 

56, -20,22 

%6  7 

58,-20,22 

%68 

49,20,20 

%5  9 

49,20,12 

%7C 

51,20,12 

%7i 

51,20,20 

%72 

52,20,20 

%73 

52,20,12 

%74 

56,20,22 

%75 

58,20,22 

%76 

60, -13,22 

%77 

60,-7,22 

%78 

60,7,22 

S-  r?  c 

60,13,22 

%8Q 

60, -7,20 

%S  1 

57,-13,20 

%S2 

60,7,20 

%83 

57,13,20 

%84 

43,-13,20 

%85 

43,13,20 

%S  6 

59,-16,20 

%87 

59, -16,10 

%88 

59,16,20 

%8S 

59,16,10 

%9G 

-39,-5,20 

%51 

-37,-4,25 

%92 

-37,4,25 

%S3 

-39,5,20 

%94 

-29,-4,25 

%95 

-28,-5,20 

-29,4,25 

497 

-28,5,20 

%9S 

-16,-3,24 

%99 

-10,-3,24 

%1C 

-10,-3,20 

%  1 

-16,-3,20 

%1C2 

-16,3,24 

%103 

-10,3,24 

%10 

-10,3,20 

Q-  ■!  n 
o  i. 

-16,3,20 

%106 
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18,-1,48 

%107 

20,-1,48 

%10S 

20,1,48 

%109 

18,1,48 

%iio 

18,-1,20 

%iii 

20, -1,20 

%112 

20,1,20 

'C  J-  J.  J. 

18,1,20 

%114 

36,-2,42 

%il5 

39,-2,42 

%116 

39,2,42 

%1 1 7 

36,2,42 

%11S 

36,-2,20 

%11 9 

39, -2,20 

%120 

39,2,20 

%121 

36,2,20 

%122 

48,-2,33 

%123 

52, -2,33 

%124 

52,2,33 

%125 

48,2,33 

%  J.  z  6 

48, -2,25 

0^  -!  O  *7 

52,-2,25 

%i28 

52,2,25 

%129 

48,2,25] ; 

%!;: 

%riumber  =  66  Dianes 


planel  =  [1, 2, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,4,4] ; 

plane2  =  [2, 5, 6, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,3,3]  ; 

plane3  =  [5, 7, 8, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6,6,6] ; 

plane4  =  [7,9,10,11,12,13,8,8,8,8,8,8,8,8,8,8,8,8]  ; 

planes  =  [12,14,15,13,13,13,13,13,13,13,13,13,13,13,13,13,13,13] ; 

planes  =  [14 , 16 , 17, 15, 15, 15 , 15 , 15, 15 , 15, 15 , 15 , 15, 15 , 15 , 15 , 15 , 15]  ,- 

plane7  =  [16,18,19,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17] ; 

plane8  =  [18,20,21,22,23,24,19,19,19,19,19,19,19,19,19,19,19,19]  ; 

plane9  =  [23,25,26,24,24,24,24,24,24,24,24,24,24,24,24,24,24,24] ; 

planelO  =  [25,27,28,26,26,26,26,26,26,26,26,26,26,26,26,26,26,26] 

planell  =  [27,1,4,28,28,28,28,28,28,28,28,28,28,28,28,28,28,28]  ; 

botton  =  [4,3,6,8,13,15,17,19,24,26,28,28,28,28,28,28,28,28] ; 

deckl  =  [1,2,27,27,27,27,27,27,27,27,27,27,27,27,27,27,27,27]; 

deck2  =  [2,5,25,27,27,27,27,27,27,27,27,27,27,27,27,27,27,27]; 

deck3  =  [5,7,9,29,30,22,23,25,25,25,25,25,25,25,25,25,25,25]; 

deck4  =  [9,10,11,31,32,33,88,88,88,88,88,88,88,88,88,88,88,88]  ; 

decks  =  [22,21,20,34,35,36,90,90,90,90,90,90,90,90,90,90,90,90]  ; 

deck6  =  [37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54] ; 

deck7  =  [12,14,16,18,18,18,18,18,18,18,18,18,18,18,18,18,18,18]  ; 

deck8  =  [55,56,57,58,59,60,60,60,60,60,60,60,60,60,60,60,60,60] ; 

laterall  =  [9,10,38,61,62,63,64,65,66,67,68,68,68,68,68,68,68,68]  ; 

lateral2  =  [22,21,51,69,70,71,72,73,74,75,76,76,76,76,76,76,76,76] 

lateral3  =  [77,78,60,59,79,80,30,29,29,29,29,29,29,29,29,29,29,29] 

lateral4  =  [68,77,29,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9]; 

laterals  =  [76,80,30,22,22,22,22,22,22,22,22,22,22,22,22,22,22,22] 

lateral6  =  [55,60,81,82,82,82,82,82,82,82,82,82,82,82,82,82,82,82] 

lateral?  =  [58,59,83,84,84,84,84,84,84,84,84,84,84,84,84,84,84,84] 

lateral8  =  [56,55,82,85,85,85,85,85,85,85,85,85,85,85,85,85,85,85] 

lateral9  =  [57,58,84,86,86,86,86,86,86,86,86,86,86,86,86,86,86,86] 

laterall0=  [56, 57, 86 , 85 , 85 , 85 , 85, 85 , 85 , 85 , 85 , 85, 85, 85 , 85 , 85 , 85 , 85]  ; 
lateralll=  [87,88,33,40,40,40,40,40,40,40,40,40,40,40,40,40,40,40]; 


laterall2= 

laterall3= 

laterall4= 

laterall5= 


[89,90,36,49,49,49,49,49,49,49,49,49,49,49,49,49,49,49] 

[40,33,32,41,41,41,41,41,41,41,41,41,41,41,41,41,41,41] 

[49,36,35,48,48,48,48,48,48,48,48,48,48,48,48,48,48,48] 

[31,42,41,32,32,32,32,32,32,32,32,32,32,32,32,32,32,32] 
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Iaterall6=  [34 , 47, 48 , 35 , 35 , 35 , 35 , 35 , 35 , 35 , 35 , 35 , 35 , 35 , 35, 35 , 35 , 35] ; 
laterall7=  [43 , 11 , 31 , 42 , 42 , 42 , 42 , 42 , 42 , 42 , 42 , 42 , 42 , 42 , 42 , 42 , 42 , 42] ; 
lateral 18=  [46,20,34,47,47,47,47,47,47,47,47,47,47,47,47,47,47,47]; 
laterall9=  [44,43,11,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12]  ; 
lateral20=  [45,46,20,18,18,18,18,18,18,18,18,18,18,18,18,18,18,18]; 
lateral21=  [44,12,18,45,45,45,45,45,45,45,45,45,45,45,45,45,45,45]; 
stackl  =  [91,92,93,94,94,94,94,94,94,94,94,94,94,94,94,94,94,94]; 

stack2  =  [91,92,95,96,96,96,96,96,96,96,96,96,96,96,96,96,96,96] ; 

stack3  =  [96,95,97,98,98,98,98,98,98,98,98,98,98,98,98,98,98,98] ; 

Stack4  =  [94,93,97,98,98,98,98,98,98,98,98,98,98,98,98,98,98,98]; 

Stack5= [92,95,97,93,93,93,93,93,93,93,93,93,93,93,93,93,93,93]; 
turretl= [99, 100, 101, 102, 102, 102, 102, 102, 102, 102, 102, 102, 102, 102, 102, 102, 102, 102 
]  ; 

turret2= [100,104,105,101,101,101,101,101,101,101,101,101,101,101,101,101,101,10 
1]  ; 

turret3= [103,104,105,106,106,106,106,106,106,106,106,106,106,106,106,106,106,10 
6]  ; 

turret4= [99, 102, 106, 103, 103, 103, 103, 103, 103, 103, 103, 103, 103, 103, 103, 103, 103, 103 

]  ; 

turret5= [99,100,104,103,103,103,103,103,103,103,103,103,103,103,103,103,103,103 

]  ; 

mastl= [107,108,112,111,111,111,111,111,111,111,111,111,111,111,111,111/111,111] 
mast2= [108,109,113,112,112,112,112,112,112,112,112,112 , 112,112,112,112,112,112] 
raast3= [109,110,114,113,113,113,113,113,113,113,113,113,113,113,113,113,113,113] 
mast4= [110 , 107,111,114,114 , 114 , 114 , 114,114,114,114 , 114 , 114 , 114,114,114 , 114 , 114] 
mast5= [107,108,109,110,110,110,110,110,110,110,110,110,110,110,110,110,110,110] 

surfacel= [115,116,120,119,119,119,119,119,119,119,119,119,119,119,119,119,119,1 

19]  ; 

surface2= [116,117,121,120,120,120,120,120,120,120,120,120,120,120,120,120,120,1 

20]  ; 

surface3= [117,118,122,121,121,121,121,121,121,121,121,121,121,121,121,121,121,1 

21]  ; 

surface4= [118,115,119,122,122,122,122,122,122,122,122,122,122,122,122,122,122,1 

22]  ; 

surface5= [115,116,117,118,118,118,118,118,118,118,118,118,118,118,118,118,118,1 
18]  ; 

radarl= [123, 124, 128, 127, 127, 127, 127, 127, 127, 127, 127, 127, 127, 127, 127, 127, 127, 127 

]  ; 

radar2= [124, 125, 129, 128, 128, 12  8, 12  8, 128, 128, 128, 128, 12  8, 128, 128, 128, 128, 128, 128 

]  ; 

radar3= [125, 126, 130, 129, 129, 129, 12  9, 12  9, 129, 129, 129, 12  9, 129, 129, 129, 129, 129, 129 

]  ; 

radar4= [126, 123, 127, 130, 130, 130, 13 0,130, 130, 130, 130, 130, 130, 130, 130, 130, 130, 130 

]  ; 

radar5= [123, 124, 125, 126, 126, 126, 126, 126, 126, 126, 126, 126, 126, 126, 126, 126, 126, 126 

]  ; 

%iaces  56  planes 
faces^  [planel 
plane2 
plane3 
plane4 
planes 
planes 
plane? 
planes 
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planes 

planelO 

planell 

botton 

deckl 

deck2 

deck3 

deck4 

decks 

decks 

deck? 

decks 

laterall 

lateral2 

laterals 

lateral4 

laterals 

laterals 

lateral? 

laterals 

laterals 

laterallO 

lateralll 

lateralis 

lateralis 

laterall4 

lateralis 

lateralis 

laterall? 

lateralis 

lateralis 

lateralSO 

lateralSl 

stackl 

Stacks 

Stacks 

stack4 

Stacks 

turretl 

turrets 

turrets 

turret4 

turrets 

mastl 

masts 

masts 

mast4 

masts 

surf acel 

surfaces 

surfaces 

surf ace4 

surfaces 

radarl 

radars 

radars 

radar 4 

radars] ; 


%  End  of  file  frigate. m 
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%  Naval  Pcstaraduace  Schcci  -  CA 


Author  :  Jorge  Aniaral  Alves,  LCDR  (Brazilian  Navy) 


function  [verts , faces] =pointsur ( )  ; 


[93, 

0,  0 

3  0 

H 

CO 

0, 

-15 

%"' 

39, 

16, 

-15 

%2 

-43 

,  16,  -15 

-43 

,“16, 

,  -15 

%4 

39, 

-16, 

-15 

%5 

51, 

-16, 

0 

%6 

47, 

-16, 

-5 

%7 

37, 

-10, 

-5 

%8 

37, 

,-16, 

12 

%9 

17, 

-16, 

12 

%10 

15, 

-16, 

-5 

%il 

-43, 

,-16, 

-8 

%12 

-43, 

16,  < 

■•8 

15, 

16, 

-5 

%14 

17, 

16, 

12 

%i5 

37, 

16, 

12 

%16 

37, 

10, 

-5 

%17 

47, 

16, 

-5 

%18 

51, 

16,  0 

%1 

9 

21, 

3, 

12 

%20 

17, 

3, 

12 

%2 1 

17, 

-3, 

12 

%22 

H 

CN 

-3, 

12 

%23 

21, 

3, 

19 

%24 

17, 

3, 

19 

%2  5 

17, 

“3, 

19 

2  6 

21, 

-3, 

19 

%27 

6, 

16, 

-8 

%28 

-6, 

16, 

-8 

%29 

-6, 

8, 

-8 

■{>  ^  0 

6, 

8, 

-8 

%3 1 

6, 

16, 

6 

%32 

-6,  16,  6  %33 

-6,  8,  6  %34 

6,  8,6  %3  5 

-16,  16,  -8  %35 

-16,  16,  -5  %37 

-16,-16,  -5  %38 

-16,-16,  -8  %3S 

38,  -16,  8  %40 

38,  16,  8  %41 

31,  16,  8  %42 

31,  16,  12  %43 
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31, 

-16, 

31, 

-16, 

17, 

-16, 

17, 

16, 

41, 

-16, 

41, 

16,  4 

33, 

16,  4 

33, 

-16, 

31, 

-16, 

31, 

16, 

31, 

16,  4 

27, 

16,  4 

25, 

16, 

29, 

16, 

25, 

16,  4 

21, 

16,  4 

19, 

16, 

23, 

16, 

31, 

-16, 

27, 

-16, 

25, 

-16, 

29, 

-16, 

25, 

-16, 

21, 

-16, 

19, 

-16, 

23, 

-16, 

37, 

10,  4 

37, 

-10, 

17, 

-10, 

17, 

-10, 

17, 

10, 

12  %44 

8  %45 

-8  %46 

-8  %47 

4  %4S 

%4  9 

%5  G 

4  %5 1 

-5  %52 

-5  %53 

%55 

-5  %56 

-5  %57 

%5S 
%59 

-5  %60 

-5  %61 

4  %62 

4  %63 

-5  %64 

-5  %65 

4  %66 

4  %67 

-  5  %  c  8 

-5  %69 

%70 

4  %71 

4  %72 

-5  %73 

-5  %74 


17,  10,  4]  ;  %75 


%nuirLber  =  32 


planeO  = 
planel  = 
plane2  = 
planes  = 
lateO  = 
latel  = 
late2  = 
lateS  = 
late4  = 
bottom  = 
glassO  = 
glassl  = 
glass2  = 
radarO  = 
radarl  = 
radar 2  = 
stackO  = 
stackl  = 
Stacks  = 
Stacks  = 
deckO  = 
deckl  = 
decks  = 
decks  = 
deck4  = 


[0,1,2,18,19,19,19,19,19,19,19,19,19,19,19,19]; 

[2.3.13.47.14.18.18.18.18.18.18.18.18.18.18.18] / 
[4,5,7,11,46,12,12,12,12,12,12,12,12,12,12,12]; 
[0,6,7,5,1,1,1,1,1,1,1,1,1,1,1,1]; 

[47,14,37,36,36,36,36,36,36,36,36,36,36,36,36,36] 
[8,71,72,73,73,73,73,73,73,73,73,73,73,73,73,73]; 
[17,70,75,74,74,74,74,74,74,74,74,74,74,74,74,74] 
[52,51,48,40,45,44, 10,11,68,67,66,69, 64,63,62,65] 
[53,50,49,41,42,43,15,14,60,59,58,61,56,55,54,53] 
[1,2, 3, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5]; 

[48,40,41,49,49,49,49,49,49,49,49,49,49,49,49,49] 
[45,44,43,42,42,42,42,42,42,42,42,42,42,42,42,42] 
[8,17,70,71,71,71,71,71,71,71,71,71,71,71,71,71] ; 

[20.24.27.23.23.23.23.23.23.23.23.23.23.23.23.23] 

[20.21.25.24.24.24.24.24.24.24.24.24.24.24.24.24] 

[22.23.27.26.26.26.26.26.26.26.26.26.26.26.26.26] 
[28,32,35,31,31,31,31,31,31,31,31,31,31,31,31,31] 

[32.33.29.28.28.28.28.28.28.28.28.28.28.28.28.28] 
[30,31,35,34,34,34,34,34,34,34,34,34,34,34,34,34] 

[30.34.33.29.29.29.29.29.29.29.29.29.29.29.29.29] 
[0,19,6,6,6,6,6,6,6,6,6,6,6,6,6,6] ; 

[6.7.18.19.19.19.19.19. 19.19.19.19.19.19.19.19]  ; 
[7,8,17,18,18,18,18,18,18,18,18,18,18,18,18,18]; 
[15,16,9,10,22,23,20,21,21,21,21,21,21,21,21,21]; 

[24.25.26.27.27.27.27.27.27.27.27.27.27.27.27.27] 
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decks  = 
decks  = 
deck?  = 
decks  = 
deck9  = 
decklO  = 
deckll  = 


[14,15,21,25,26,22,10,11,11,11,11,11,11,11,11,11] 

[11,38,37,29,30,31,28,14,14,14,14,14,14,14,14,14] 

[36.37.38.39.39.39.39.39.39.39.39.39.39.39.39.39] 

[36.13.12.39.39.39.39.39.39.39.39.39.39.39.39.39] 
[3,4,12,13,13,13,13,13,13,13,13,13,13,13,13,13]/ 
[32,33,34,35,35,35,35,35,35,35,35,35,35,35,35,35]; 
[40,45,42,41,41,41,41,41,41,41,41,41,41,41,41,41] ; 


faces= [  planeO 
planel 
plane2 
plane3 
lateO 
latel 
late2 
late3 
late4 
bottom 
glassO 
glassl 
glass2 
radarO 
radarl 
radar2 
stackO 
stackl 
stack2 
Stacks 
deckO 
deckl 
deck2 
decks 
deck4 
decks 
decks 
deck? 
decks 
deck9 
decklO 
deckll] ; 


faces=faces+l  ,* 


%  End  of  file  point sur.m 
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Naval  Postgraduate  School  -  CA 


Bate  01 
Version 
Author  : 


Function 
merchant . m 

:  returns  the  vertices  and  faces  to  be  used  by  the  iyATLA.3 
’■patch”  in  order  to  construct  the  merchant  3-D  model 
march  2001 
:  1.0 

Jorge  Asmara  1  Alves,  LCDR  (Brazilian  Navy) 


:ion 


function  [verts , faces] =merchant ( ) ; 


verts=  [115 , -11 , 0 

%1 

115,-15,10 

%2 

118.9,-14.5,10.8 

%3 

117.8,-10.6,0 

%4 

122.5,-13,11.5 

fee 

120.5,-9.5,0 

%6 

125.6,-10.6,12 

%7 

122.8,-7.8,0 

%8 

128,-7.5,12.5 

%9 

124.5,-5.5,0 

%1C 

129.5,-3.9,13 

%  X 

125.6,-2.85,0 

%12 

130,0,13 

%13 

126,0,0 

%14 

129.5,3.9,13 

%15 

125.6,2.85,0 

%16 

128,7.5,12.5 

%17 

124.5,5.5,0 

%18 

125.6,10.6,12 

%  19 

122.8,7.8,0 

%2  0 

122.5,13,11.5 

%21 

120.5,9.5,0 

%22 

118.9,14.5,10.8 

%23 

117.8,10.6,0 

%24 

115,15,10 

%25 

115,11,0 

%2  6 

96,-15,10 

%27 

93,-15,7 

%28 

-130,-15,7 

%29 

-128,-11,0 

%'3  0 

-128,11,0 

%31 

-130,15,7 

%32 

93,15,7 

%33 

96,15,10 

%34 

-80,-15,23 

%3  5 

-80,15,23 

%3  6 

-84,15,23 

%37 

-84,11,23 

%3S 

-91,11,23 

%39 

-91,-11,23 

%40 

-84,-11,23 

%41 

-84,-15,23 

%42 

-91,-11,17 

%43 

-91,11,17 

%44 

-110,11,17 

%45 

-110,-11,17 

%46 
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-110,11, 7 

^47 

-110, -11,7 

%4 

-80,-15,21 

%49 

-80,-11,19 

%50 

-80,-11,7 

% 1 

-80,11,7 

%  5  2 

-80, 11,19 

%53 

-80,15,21 

%54 

-84, -15,21 

%55 

-84,-11, 19 

%56 

-91, -11,7 

%  5  7 

-84,15,21 

%5  S 

-84,11,19 

%5  9 

-91,11,7 

^60 

-99,-3,17 

%61 

-99,3,17 

%62 

-104,3,17 

%  6  3 

-104,-3,17 

%64 

-99,3,27 

%65 

-99, -3,27 

6 

-104,3,27 

%6  7 

-104, -3,27 

%6S 

2,  -11,7 

%6  9 

2, -9, 7 

%7  0 

0,-9,7 

0,-11, 7 

%72 

2,-11,20 

%73 

2, -9,20 

%74 

0,-11,20 

%75 

0, -9,20 

%76 

2,9,7 

%77 

2,11,7 

%78 

0,11,7 

%79 

0,9,7 

%80 

2,9,20 

%S1 

2,11,20 

%S2 

0,11,20 

%S3 

0,9,20 

^84 

111,-1,10 

%85 

111,1,10 

%36 

109,1,10 

'<5  87 

109, -1,10 

%88 

111,-1,26 

%  8  9 

111,1,26 

%90 

109,1,26 

^  Q  -r 

109,-1,26 

%92 

-89, -1,23 

%S3 

-89,1,23 

%  Q4 

-91,1,23 

%95 

-91, -1,23 

%  9  6 

-89, -1,30 

%97 

-89,1,30 

%9S 

-91,1,30 

%99 

-91,-1,30]  ; 

%  1 0( 

%Dlanes 


proal  =  [1, 2, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,4,4], • 
proa2  =  [3,5,6,4,4,4,4,4,4,4,4,4,4,4,41; 
proa3  =  [5,7,8,6,6,6,6,6,6,6,6,6,6,6,61; 
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proa4  =  [7,9,10,8,8,8,8,8,8,8,8,8,8,8,8]; 

proas  =  [9,11,12,10,10,10,10,10,10,10,10,10,10,10,10]; 

proa6  =  [11,13,14,12,12,12,12,12,12,12,12,12,12,12,12]; 

proaV  =  [13,15,16,14,14,14,14,14,14,14,14,14,14,14,14]; 

proas  =  [15,17,18,16,16,16,16,16,16,16,16,16,16,16,16]; 

proa9  =  [17,19,20,18,18,18,18,18,18,18,18,18,18,18,18]; 

proal0=  [19,21,22,20,20,20,20,20,20,20,20,20,20,20,20] ; 
proall=  [21,23,24,22,22,22,22,22,22,22,22,22,22,22,22] ; 
proal2=  [23,25,26,24,24,24,24,24,24,24,24,24,24,24,24] ; 
planel  =  [1,2,27,28,29,30,30,30,30,30,30,30,30,30,30]; 

plane2  =  [29,30,31,32,32,32,32,32,32,32,32,32,32,32,32]; 

plane3  =  [26,25,34,33,32,31,31,31,31,31,31,31,31,31,31]; 

deckl  =  [2,3,5,7,9,11,13,15,17,19,21,23,25,25,25]; 

deck2  =  [2,25,34,27,27,27,27,27,27,27,27,27,27,27,27]; 

decks  =  [27,34,33,28,28,28,28,28,28,28,28,28,28,28,28]; 

deck4  =  [28,33,32,29,29,29,29,29,29,29,29,29,29,29,29]; 

decks  =  [35, 36, 37, 38, 39, 40, 41, 42, 42, 42,-42, 42, 42, 42, 42]  ; 
deck6  =  [40,39,44,43,43,43,43,43,43,43,43,43,43,43,43]; 

deck7  =  [43,44,45,46,46,46,46,46,46,46,46,46,46,46,46]; 

decks  =  [46,45,47,48,48,48,48,48,48,48,48,48,48,48,48]; 

botton  =  [1,4,6,8,10,12,14,16,18,20,22,24,26,31,30]; 

front  =  [35,49,50,51,52,53,54,36,36,36,36,36,36,36,36]; 

laterall  =  [51,50,56,41,40,43,57,57,57,57,57,57,57,57,57]; 

lateral2  =  [52,53,59,38,39,44,60,60,60,60,60,60,60,60,60]; 

laterals  =  [43,46,48,57,57,57,57,57,57,57,57,57,57,57,57]; 

lateral4  =  [44,45,47,60,60,60,60,60,60,60,60,60,60,60,60]; 

blocorightl  =  [35,49,55,42,42,42,42,42,42,42,42,42,42,42,42]; 

blocoright2  =  [41,42,55,56,56,56,56,56,56,56,56,56,56,56,56]; 

blocoleftl  =  [36,37,58,54,54,54,54,54,54,54,54,54,54,54,54]; 

blocoleft2  =  [37,38,59,58,54,54,54,54,54,54,54,54,54,54,54]; 

stackl  =  [64,61,66,68,68,68,68,68,68,68,68,68,68,68,68]; 

Stacks  =  [61,62,65,66,66,66,66,66,66,66,66,66,66,66,66]; 

Stacks  =  [62,65,67,63,63,63,63,63,63,63,63,63,63,63,63]; 

stack4  =  [67,68,64,63,63,63,63,63,63,63,63,63,63,63,63]; 

stacktop  =  [66,65,67,68,68,68,68,68,68,68,68,68,68,68,68]; 

middlerightl  =  [72,69,73,75,75,75,75,75,75,75,75,75,75,75,75]; 

middleright2  =  [69,70,74,73,73,73,73,73,73,73,73,73,73,73,73]; 
middleright3  =  [70,74,76,71,71,71,71,71,71,71,71,71,71,71,71]; 

middleright4  =  [71,76,75,72,72,72,72,72,72,72,72,72,72,72,72]; 

middlerighttop  =  [73,74,75,76,76,76,76,76,76,76,76,76,76,76,76]  ; 

middleleftl  =  [80,77,81,84,84,84,84,84,84,84,84,84,84,84,84]; 

middleleft2  =  [77,78,82,81,81,81,81,81,81,81,81,81,81,81,81]; 

middlelefts  =  [78,79,83,82,82,82,82,82,82,82,82,82,82,82,82]; 

middleleft4  =  [79,80,84,83,83,83,83,83,83,83,83,83,83,83,83]; 

middlelefttop  =  [81,82,83,84,84,84,84,84,84,84,84,84,84,84,84]  ; 

mastproal  =  [88,85,89,92,92,92,92,92,92,92,92,92,92,92,92]; 

mastproaS  =  [85,86,90,89,89,89,89,89,89,89,89,89,89,89,89]; 

mastproaS  =  [86,87,91,90,90,90,90,90,90,90,90,90,90,90,90]; 

mastproa4  =  [87,88,92,91,91,91,91,91,91,91,91,91,91,91,91]; 

mastproatop  =  [89,90,91,92,92,92,92,92,92,92,92,92,92,92,92]; 

mastpopal  =  [96,93,97,100,100,100,100,100,100,100,100,100,100,100,100]  ; 

mastpopas  =  [93,94,98,97,97,97,97,97,97,97,97,97,97,97,97]; 

mastpopaS  =  [94,98,99,95,95,95,95,95,95,95,95,95,95,95,95],- 

mastpopa4  =  [95,99,100,96,96,96,96,96,96,96,96,96,96,96,96]; 
mastpopatop  =  [97,98,99,100,100,100,100,100,100,100,100,100,100,100,100] 

%f aces 

faces= [proal 
proa2 
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proa3 

proa4 

proaS 

proa6 

proa? 

proaS 

proa9 

proalO 

proall 

proal2 

planel 

plane2 

planes 

deckl 

deck2 

decks 

deck4 

decks 

deck6 

deck? 

decks 

botton 

front 

laterall 

laterals 

laterals 

lateral4 

blocorightl 

blocorightS 

blocolef tl 

blocolef t2 

stackl 

Stacks 

Stacks 

stack4 

stacktop 

middlerightl 

middlerightS 

middlerightS 

middleright4 

middle right top 

middlelef tl 

middlelef t2 

middlelef tS 

middlelef t4 

middlelef ttop 

mastproal 

mastproaS 

mastproaS 

mastproa4 

mastproatop 

mastpopal 

mastpopaS 

mastpopaS 

mastpopa4 

mastpopatop] ; 


End  of  file  merchant. m 
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Naval  ir^c  St  graduate  School  -  CA 

Type  :  Function 

Name  :  iindInputMomSet .m 

Function  :  Finds  the  moment  invariants  input  set  for  each  view  angle 

This  set  v/ill  be  used  as  the  input  for  the  neural  netv/ork  classifier 
Date  01  march  2001 
Version  ;  1.0 

Author  :  Jorge  Amaral  Alves,  LCDR  (Brazilian  Na\/"y) 


function  [trainingSet]  =  f indInputMomSet (  noise_silhouettes, rows, columns) ; 

totalLines=12 ;  %number  of  moment  functions 
[X, totalColumns] =si2e (noise_silhouettes) ; 
trainingSet=  zeros (totalLines , totalColumns) ; 

Xl=zeros (rows, columns) ; 

for  i=l : totalColumns 

temp=noise_silhouettes ( : , i) ; 

XI ( : ) =temp; 

trainingSet (1 :6, i) =find_mom_f unctions (XI) ; 

X2=X1*255/ 

XXl=edge (X2 , ' prewitt ‘ ) ; 

0  ...  1  max-1 

trainingSet (7 : 12 , i) -find_mom_f unctions (XXI) ; 

end 


%  End  of  file  f indInputMomSet .m 


%  solid  silhouette 

% [  0  0  . . , i  max=255 
%only  the  edges  [0  0 

%boundary 


% 

% 
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Maval  Postgraduate  School  -  CA 

Type  :  Function 

Name  :  f  ind__mom__f unction s  .m 

Function  :  returns  six  functions  values  relating  to  the  central  moment 
notions 

invariant  under  rotation,  translation, ref lection  and  scale 
Date  01  march  20C1 
Version  :  1.0 

Author  :  Jorge  Amaral  Alves,  LCDR  (Brazilian  Navy) 


function  NewMoments  =  find_mom_functions (pattern) ; 

IND_MOM_FUNCTIONS  returns  six  functions  values  relating  to  the  central  moment 
notions 

nvariant  under  rotation,  translation, ref lection  and  scale 

[Nsw.^2  ,  Ne'A'M3  ,  Nev;M4  ,  KewM5  ,  NewM6 ,  Nswl^7 ]  =  f  ind_mom_f unctions  (pattern) 

%find  the  second  and  third-order  central  moments 
ind_moment (1, 1, pat tern) ; 
m_l_2=f ind_moment ( 1 , 2 , pattern) ; 
m_2_l=f ind_moment (2 , 1 , pattern) ; 
m_2_0=f ind_moment (2 , 0 , pattern) ; 

Tn_0_2=f  ind_moment  ( 0 , 2 ,  pattern)  ; 
m_0_3=f ind_moment (0,3, pattern) ; 
m_3_0=f ind_moment (3, 0, pattern) ; 

%iind  the  moment  functions  invariant  under  rotation  and  reflection 
M2=(m_2_0  -  m_0_2)^2  +  4*(m_l_1^2); 

M3=(m_3_0  -  3*m_l_2)"2  +  (3*m_2_l  -  m_0_3)^2; 

M4=(m_3_0  +  m_l_2)^2  +  (m_2_l  +m_0_3)^2; 

M5=(m_3_0  -  3*m_l_2) * (m_3_0  +  m_l_2 ) * ( {m_3_0  +  m_l_2)"2  -  3*(m_2_l  + 

m_0_3) "2)  +  (3*m_2_l  -  m_0_3 ) * (m_2_l  +  m_0_3 ) * ( 3  * ( (m_3_0  +  m_l_2)"2)  -  (m_2_l  + 

m_0_3) ^2) ; 

M6=(m_2_0  -  m_0_2)  *  (  (m_3_0  +  tn_l_2)^2  -  (Tn_2_l  +  m_0_3)"2)  +  4*m_l_l*  (m_3_0  + 

m_l_2) * (m_2_l  +  m_0_3 ) ; 

M7=(3*m_2_l  -  in_0_3)  *  (m_3_0  +  m_l__2 )  *  (  (m_3_0  +  m_l_2)^2  -  3*(m_2_l  +  m_0_3)^2) 
-  (m_3_0  ~  3*m_l_2) * (m_2_l  +  m_0_3 ) * (3* (m_3_0  +  m_l_2)"2  -  (m_2_l  +  m_0_3)^2); 

%normalizing  the  moment  functions  under  scale  using  the  radius  of  gyration 
r=sqrt  (m_2_0  +  i^_0_2)  ; 

%the  new  moment  function  M2  through  M7 
NewM2  =  M2/ {r"4) ; 

NewM3  =  M3/ (r^6) ; 

NewM4  =  M4/ (r^7) ; 

NewMS  M5/(r^l2); 

NewM6  =  M6/ (r"8) ; 

NewM7  =  M7/ (r^l2) ; 

NewMoments= [NewM2 
NewM3 
NewM4 
NewMS 
NewM6 
NewM7] ; 


%  End  of  file  find  mom  functions. m 
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Naval  Pcsccrraduate  School  - 


%  Type  :  Function 

%  Name  :  f  ind_moTnent .  m 

%  Function  :  returns  the  central  moment  related  to 
%  Date  01  march  2001 

%  Version  :  1.0 

1  Author  :  Jorge  Amaral  Alves,  LCDR  (Brazilian  Navy) 


the  indexes  p  and  q 


function  [in_p_q]  =  find_inoment  (p, q, pattern)  ; 

%FIND_MOMENT  returns  the  central  moment  related  to  the  indexes  p  and  q. 
%  =  ?IND_MOMENT  (  ? ,  Q ,  PATTERN) 

[NRows, NColumns] =size (pattern) ; 
count = sum  (sum (pat tern,  1)  ,2) 

'5 rind  tne  centroid  values  nor  tnis  pattern 


[um,vTn]=  find_centroid  (pattern)  ; 

%index  rov;s==>vi  index  columns==>ui 
^  u  i  -  um )  ^  p }  *  (  ( V i  -  vm )  ^  q ) 

[ui ,  vi]  =:meshgrid  ( 0  :  (NColumns-1)  ,  0  :  (NRows-1)  )  ;  %ui  = 

% 


% 


[0 

n 


0 


2 

2 


2  2  2  .  .  . 


m_p_q= ( (ui-um)  . ^p)  . * ( (vi~vm)  . ^q)  ; 
pattern=double (pattern) ; 
m_p_q=m_p_q . *pat  t  ern ; 
m_p_q=  sum { sum (m_p_q , 1 ) , 2 ) ; 
m_j)_q=m_p_q/  count  ; 

—  —  —  —  —  —  —  —  —  =  =  = 
%  End  of  file  finf  rnoinent.m 
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%  Naval  Postgraduate  School  -  CA 


%  Type  :  Function 

%  Name  :  find  centroid.m 


Author  :  Jorge  Amaral  Alves,  LCDR  (Brazilian  Navy) 


function  [um,vm]  =  find_centroid (pattern) ; 

%?IND_CENTROID  returns  the  centroid  values  for  an  image  frame  PATTERN. 
%  [VM,UMj  =  ?IND_CENTROID( PATTERN) 


%considering  the  solid  silhouette  or  its  boundary  as  a  binary  matrix (ones  and 
zeros) 

pattern=double (pattern) ; 

[NRows^NColumns] = size (pattern) ; 

%th6  mean  values  um  and  vm  are  the  centroid  of  the  civen  oattern 


vi=  [0  0  C  .  .  . 
111... 
2  2  2  .  .  . 


%  index  rov7S  =  =  >vi  index  cclumns==:>ui 

(  (ui-um)  ^p)  *  (  (vi-umn)  ^q) 

[ui,vi] =meshgrid(0 : (NColumns-1) , 0 : (NRows-1) ) ;%ui= [0  i  2 

%  0  12 

%  0  12 


ui=ui . *pattern; 

um=  sum ( sum ( ui , 1 ) ,2)  ; 

vi=vi . *pattern; 

vm=sum ( sum ( vi , 1 ) ,2)  ; 

count = sum ( sum (pattern, 1) ,2) ; 

vm=vm/ count ; 

um=um/ count ; 

%  End  of  file  find  centroid.m 
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%  Naval  Postgraduate  School  - 

%  Type  :  Main  program 

%  Name  :  mainShipRecon 

%  Function  :  Create  and  trains  a  neural  n8tv;ork  responsible  for  recognizing 

%  based  on  the  moment  invaz'iants  calculated  for  each  viewangle 

silhouette 

%  Date  Cl  march  2001 

%  Version  ;  i .  0 

%  Author  :  Jorge  Amaral  Alves,  LCDR  (Brazilian  Navy) 


azAngleInc=15 ;  %  ==>12  azimuth  angles 

eleAngleInc=15 / %0 , 15 , 3C , 45  ==>four  elev  angles  %  9C;%one  ele\“aticn  angle  ==> 
zero  degree 

totalAzimuth=180/azAngleInc; %  ==>12 
totalElevation=4 ; %=90/eleAngleInc; 

totalColumns=totalAzimuth*totalElevation;%each  column  means  one  viev; 
totalLines=12 ;  %number  of  moment  functions 

inputSet=  zeros (totalLines, totalColumns) ; 
rows=210 ; 
coluinns=2  80  ; 

silhouettes=zeros (rows* columns , total Columns) ; 


DEFINING  THE  MODEL  PROBLEM 


The  script  file  FINDINPUTMOMSET  defines  a  matrix  inputSet 
which  contains  the  moment  functions (12  \^alues)  for  all  views  defined  by 

"Azimuth  Angle  Increment"  and  "Elevation  .Angle  Increment"  of  the  ship 


Each  target  vector  has  5  elements  Vvdth 

all  zeros,  except  for  a  single  i.  Aircraft  Carrier  has  a  1  in  the 
first  element,  Destroyer  in  the  second,  Frigate  in  the  third 
,  Point  Sur  in  the  fourth  and  Merchant  in  the  fifth. 


%aircraft  carrier 

[verts , faces] =aircarrier; 

[inputSet , silhouettes] =  f indInputSet (  verts, faces, azAngleInc,eleAngleInc) ; 
[x, y] =size ( inputSet) ; 

temp=zeros {5,y) ; 

temp (1, : ) =1 ; %f irst  line  means  aircraft  carrier 


% 


% 

% 

the 

model 


targets=temp; 
alphabets inputSet ; 
all_silhouettes=silhouettes ; 

1 

% destroyer 

[verts , faces] =destroyer ; 

[inputSet , silhouettes] =  f indInputSet (  verts, faces, azAngleInc, eleAngleInc) ; 
[x, y] =size { inputSet ) ; 
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temp=zeros ( 5 , y) ; 

temp {2 ,:) =1 ; %5eccnd  line  means  destroyer 

targets= [targets, temp] ; 
alphabet^ [alphabet , inputSet]  ; 

all_silhouettes= [all_silhouettes, silhouettes]  ; 

2 

Irrigate 

[verts , faces] =f rigate ; 

[inputSet, silhouettes] =  f indInputSet (  verts, faces, azAngleInc, eleAngleInc) 
[x,y] =size (inputSet) ; 
temp^zeros (5 , y) ; 

temp (3 , : ) =1 ; Ithird  line  means  frigate 

targets= [targets , temp] ; 
alphabet= [alphabet , inputSet]  ; 

all_silhouettes= [all_silhouettes , silhouettes]  ; 

3 

%point  sur 

[verts , faces] =pointsur ; 

[inputSet, silhouettes] =  f indInputSet (  verts, faces , azAngleInc , eleAngleInc) 
[x, y] =size ( inputSet ) ; 
temp=zeros ( 5 , y) ; 

temp (4 , : ) =1 ; %f ourth  line  means  point  sur 

target s=  [targets ,  temp] 
alphabet= [alphabet, inputSet] ; 

all_silhouettes= [all_silhouettes , silhouettes]  ; 

4 

%merchant 

[verts , faces] ^merchant ; 

[inputSet , silhouettes] =  f indInputSet (  verts , faces , azAngleInc , eleAngleInc) 
[x, y] =size (inputSet) ; 
temp=zeros ( 5 , y) ; 

temp  ( 5  ,  : )  =1 ;  %f  if  th  line  m.eans  point  sur 

targets^ [targets , temp] ; 
alphabet= [alphabet , inputSet] ; 

all_silhouettes= [all__silhouettes , silhouettes]  ; 

save  all_silhouettes  all_silhouettes / 
clear  all_silhouettes ; 
save  alphabet  alphabet  ,- 
save  targets  targets; 

Clf  ; 

figure (gcf ) 
echo  on 

[R,Q]  =  size (alphabet) ; 

[S2,Q]  =  size (targets) ; 

%pause  %  Strike  any  key  to  define  the  network... 


%  DEFINING  THE  NETWORK 
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The  ship  recognition  netVv^crk  v/ili  have  20  TANSIG 
neurons  in  its  hidden  layer. 


SI  =  20; 

net  =  newff (minmax (alphabet ), [SI  S2] , { ’ logsig '  ' logsig ' } , ’ traingdx' ) ; 

net.LW{2,l}  =  net.LW{2,l}*0.01; 
net.b{2}  =  net .b{2 } *0 . 01 ; 

%pause  %  Strike  any  key  to  train  the  netv;ork... 


TRAINING  THE  NETWORPC 


net .perf ormFcn  =  'sse';  %  Sum-Squared  Error  performance  function 

net . trainParam. goal  =0.1;  %  Sum -squared  error  goal. 

net .trainParam. show  =  10000 ;%20;  %  Frequency  of  progress  displays 

epochs) , 

net . trainParam. epochs  =  600000/  %  Maximum  number  of  epochs  to  train, 

net . trainParam. me  =0.95;  %  Momentum  constant. 

%  Training  begins.  .  .please  -wait.  .  . 

P  =  alphabet; 

T  =  targets; 

[net,tr]  =  train (net , P, T) ; 
save  net  net 
5 

%  ...and  finally  finishes. 

%  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

%  End  of  file  mainShipP.econ.m 
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Naval  Postgraduate  School 
2.^'- 


Name  : 
Functio: 


1  ;  Creates  the  5-D  model,  change  the  ^ 

extract  the  silhouette  for  each  trai 
Function 
0 / mav/ 2  0  0  0 


Le^vpoint  of  the  3-D  model, 


Version  :  l.C 
Author  :  Joroe  Amaral 


Ives,  LCDR  (Brazilian  Navy; 


function  [trainingSet , silhouettes]  =  f indInputSet ( 

verts, faces, azAngleInc, eleAngleInc) ; 

totalAzimuth^lSO/azAnglelnc ; 
totalElevation=4 ; %9G /eleAngleInc ; 

totalColumns=:totalAzimuth*totalElevation; leach  column  means  one  viev; 

totalLines^l2 ;  Inumber  of  moment  functions 

trainingSet=  zeros (totalLines, totalColumns) ; 

rows=210; 

columns=280 ; 

silhouettes=zeros ( rows* col umn s , totalColumns) ; 
figNumber= figure (  ... 

' Name Silhouette ’  ,  ... 

'  Position [120  120  280  210]);%  matriz  =>  210rov;s  Z  280columns 


figure (f igNumber) 

ICreates  the  3-D  model 

patch ( ’Vertices ’ , verts, ' Faces ' , faces) ; 

view (3 ) 

axis  equal ; 

axis  off 

axis  vis3d 

lineCount=0 ; 


for  i=-90 : azAngleInc : 89 
i 

%for  j  =  0  :  eleAjtgleInc  :  39 
for  j =0 : eleAngleInc : 46 

j 

lineCount=:lineCount  +  l  ; 


%azimuth  popa  until  proa 

%elevation 

%elevation 


%Changes  the  model  viewpoint  to  the  desired  training 
view{i, j ) ; 

[X, map] = capture;  %[65,  65,  ...] 

Xl=X-65;  %  [0  0  ...  j  max=l 

silhouettes ( : , lineCount ) =X1 ( : ) ; 

trainingSet (1:6, lineCount) =find_mom_f unctions (XI) ; 
X2=X1*255; 

XXl=edge (X2, ’prewitt ' ) ; 

0  , . . ]  max=l 

trainingSet (7:12, lineCount) =find_mom_f unctions (XXI) ; 

end 

end 


angl 


e 


%  solid  silhouette 
% [  0  0  . . . ]  max=255 
%only  the  edges  [0  0 

Iboundary 


-0:=  =  =  =:=  =  =  =  =  =:  =  =::  =  =  —  =  =  =  =  =  =  =  =:=:  =  =  =  : 

%  End  of  file  f indInputSet .m 
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function  interface (action) 

if  nargin<l, 

action= ' initialize ' ; 

end; 

if  strcmp (action, ' initialize ') , 
oldFigNuniber=watchon  ; 

figNumber= figure (  ... 

'Name 'Neural  Network  Ship  Silhouette  Recognition', 
'NumberTitle ' , 'off ' ,  ... 

' Visible ' , ' of f ' ,  ... 

' Backings tore ' , ' of f ' ) ; 

axPos=[0.40  0.95-0.28  0.20  0.28]/ 
axHndll=axes (  ... 

' Units ' , ' normalized ' ,  ... 

' Position axPos ,  ... 

'XTick'  ,  []  ,  '  YTick'  ,  []  ,  ... 

' Box ' , ' on ' ) ; 

labelPos= [0.05  0.80  0.30  0.05]; 
uicontrol (  ... 

' Style ' , ' text ' ,  ... 

' String ', 'Original  Ship  Silhouette',  ... 

' BackgroundColor '  ,  'k' ,  ... 

' ForegroundColor ' , 'w' ,  ... 

'Units normalized ' ,  ... 

' Position ' , labelPos) ; 

axHndl2=axes (  ... 

' Uni ts normalized ' ,  ... 

' Position axPos+ [ 0  -0.31  00],  ... 

'XTick'  ,  []  ,  'YTick'  ,  []  ,  ... 

' Box ' , ' on ' ) ; 

labelPos= [0.05  0.50  0.30  0.05]; 
uicontrol (  ... 

'Style' , 'text' ,  ... 

' String Ship  Silhouette  with  Noise',  ... 
'BackgroundColor' , 'k' ,  ... 

' ForegroundColor ' , ' w ' ,  ... 

'Units normalized' ,  ... 

' Position ' , labelPos) ; 
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axHndl3=axes (  ... 

'Units’ , 'normalized' ,  ... 

' Position axPos+ [0  -0.62  0  0] , 
'XTick'  ,  []  ,  'YTick'  A]  ,  -  .  . 

' Box ' , ' on ' ) ; 

labelPos= [0.05  0.20  0.30  0.05]; 
uicontrol (  ... 

’ Style ' , ' text ' ,  ... 

'  String Netv/ork ’’ s  Guess',  .. 
' BackgroundColor ' , 'k' ,  ... 

' ForegroundColor ' , ' w ' ,  ... 

' Units normalized ’ ,  ... 

'Position' ,labelPos) ; 


=  =  =  ^  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

%  Information  for  all  buttons 
top=0 . 95 ; 
bottom=0 . 05 ; 

labelColor= [0.8  0.8  0.8]; 

btnWid=0 . 20 ; 

btnHt=0.10; 

right=0 . 95 ; 

lef  t=right-btnWid; 

%  Spacing  betv/een  the  button  and  the  next  command's  labe 
spacing=0 . 05 ; 


%  The  CONSOLE  frame 
f rmBorder=0 . 02 ; 

frmPos= [lef t-frmBorder  bottom- frmBorder 

0 . 9+2*f rmBorder]  ; 

h=uicontrol (  ... 

' Style ' , ' frame ' ,  ... 

’  Units normalized ' ,  ... 

' Position *, frmPos,  ... 

' BackgroundCol or ' , [ 0 . 5  0.5  0.5]); 


%  The  NEW  SILHOUETTE  button 
btnNuTnber=l  ; 

yPos=top-btnHt- (btnNumber- 1) * (btnHt+spacing) ; 
labelStr= ' New  Silhouette'; 
callbackStr= ' interface { ' ' new 

%  Generic  button  information 
btnPos=[left  yPos  btnWid  btnHt] ; 
f lyHndl=uicontrol {  ... 

' Style 'pushbutton' ,  ... 

'Units' , 'normalized' ,  ... 

' Position btnPos,  ... 

' String labelStr,  ... 

'Callback' , callbackStr) ; 


%  The  NOISE  slider 

btnNumber =2 ; 

yPos=top-btnHt- (btnNumber- 1) * (btnHt+spacing) ; 
labelStr= 'Noise ' ; 

%  Generic  button  information 


btnWid+2  *  f  rmBorder 
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sldPos= [left  yPos  btnWid  btnHt/2] ; 
labelPos= [lef t  yPos+btnHt/2  btnWid  btnHt/2]; 
sldHndl^uicontrol (  ... 

' Style ' , ' slider ' ,  ... 

' Units ' , ' normalized ’ ,  ... 

’ Position' , sldPos) ; 

uicontrol (  ... 

' Style ’ , ' text ' ,  ... 

'Units ‘normalized’ ,  ... 

' String Noise ’ ,  ... 

' Position ' , labelPos)  ; 


%  The  ROTATION  slider 
btnNumber=3 ; 

yPos=top-btnHt- (btnNumber~l) * (btnHt+spacing) 
labelStr= ’Rotation ' ; 

%  Generic  button  information 

rotPos= [left  yPos  btnWid  btnHt/2]; 
labelPos=[left  yPos+btnHt/2  btnWid  btnHt/2]; 
rotHndl=uicontrol {  ... 

’ Style ’ , ' slider' ,  ... 

' Units normalized ' ,  ... 

' Position’ , rotPos) ; 

uicontrol {  ... 

’ Style ’ , ‘ text ’ ,  ... 

' Units normalized ' ,  ... 

' String ' , ' Rotation ’ ,  ... 

' Position • , labelPos)  ; 


%  =  =  =  =  =  =  =  =  =  =  =  =  ^  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

%  The  SCALE  slider 
btnNuTnber=4  ; 

yPos=top-btnHt- (btnNumber-l) * (btnHt+spacing) 
labelStr= ' Scale ' ; 

%  Generic  button  information 

sclPos=[left  yPos  btnWid  btnHt/2]; 
labelPos= [left  yPos+btnHt/2  btnWid  btnHt/2]; 
sclHndl=uicontrol (  ... 

’ Style ' , ’ slider' ,  ... 

' Units ' , ' normalized ' ,  ... 

' Position' , sclPos) ; 

uicontrol (  ... 

' Style ' , ' text ' ,  ... 

' Units normalized ' ,  ... 

' String ' , ' Scale ’ ,  ... 

'Position' , labelPos)  ; 


%  The  INFO  button 
labelStr^ ' Info ' ; 

callbackStr= ' interface ( ' ' info ' ' ) ’ ; 
inf oHndl=uicontrol (  ... 

' Style ' , ' push ’ ,  ... 
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'  Unit s normalized '  ,  ... 

'  Position [left  bottom+btnHt+spacing  btnWid  btnHt] ,  . 

’  String ' , labelStr ,  . . . 

'Callback' , callbackStr) ; 


%  The  CLOSE  button 
labelStr= ' Close ' ; 
callbackStr= ‘ close (gcf) '; 
closeHndl=uicontrol (  ... 

' Style ' , ' push ' ,  ... 

‘Units' , 'normalized' ,  ... 

'Position' , [left  bottom  btnWid  btnHt] ,  ... 

' String ' , labelStr ,  . . . 

'Callback' , callbackStr) ; 

%  Uncover  the  figure 

hndlList= [axHndll  axHndl2  axHndlS  sldHndl  rotHndl  sclHndl] 
set {f igNumber,  . . . 

'Visible ' , ' on ' ,  ... 

'UserData' ,hndlList) ; 

watchof f (oldFigNumber) ; 
choice=0 ; 
interface  new; 
f igure (f igNumber) ; 


elseif  strcmp (action, 'new ') , 
f igNumbe  r- wa  t  chon ; 

hndlList=get (f igNumber, 'Userdata') ; 
axHndll=hndlList (1); 
axHndl2=hndlList (2); 
axHndl3=hndlList (3) ; 
sldHndl=hndlList (4) ; 
rotHndl =hndlList (5)  ; 
sclHndl=hndlList (6) ; 

load  all_silhouettes; 
load  rows; 
load  columns ; 
load  net ; 

[X,  Y]  =size (all_silhouettes)  ; 
load  choice; 
choice=choice+l ; 
save  choice  choice,- 
choice 

%find  a  random  silhouette  as  the  original  one 
randsilhouette=all_silhouettes { : , choice) ; 

%firid  the  corresponding  silhouette  after  rotation 

rotsilhouette=randsilhouette ; 

rotationlevel=get (rotHndl, 'Value' ) ; 

rotationlevel*180 

Xl^zeros  (rows,  columns)  ,- 

XI ( : ) =randsilhouette; 

X2=imrotate (XI , (rotationlevel*180) , ' bilinear ' , ' crop ' ) ; 
rotsilhouette=X2 ( : ) ; 
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%find  the  corresponding  silhouette  after  scaling 

sclsilhouette=rotsilhouette; 

scalelevel=get (sclHndl , 'Value' ) 

XI = zeros (rows, columns) ; 

XI ( : ) =rotsilhouette; 

X2=imrotate (XI, 45, 'bilinear' , 'crop' ) ; 
sclsilhouette=X2 ( : ) ; 

%find  the  corresponding  silhouette  after  adding  noise 

noiselevel=get (sldHndl, 'Value') 

noise=round (randn (rows*columns, 1) *noiselevel) ; 

y=noise'>0  ,- 

w=noise<0 ; 

noise=y+w; 

noisesilhouette=rotsilhouette+noise ; 
testsilhouette=noisesilhouette ; 

testMoments=  f indInputMomSet (  testsilhouette, rows, columns) 

A  =  sim{net ,  testMoments)  ,- 

output=compet (A) ; 

result=find(output==l) -1; 

result 

outsilhouette=all_silhouettes ( : , 48*result+25) ; 
axes (axHndll) ; 

plotSilhouette ( rands ilhouette, rows, columns) ; 
axes {axHndl2) ; 

plotSilhouette (testsilhouette, rows, columns) ; 
axes (axHndl3) ; 

plotSilhouette (outsilhouette, rows, columns) ; 
watchoff (figNumber) ; 

elseif  strcmp (action, ' info ') , 
ttlStr=get (gcf , 'Name' ) ; 
hlpStr=  { ' 

'  This  window  demonstrates  the  use  of  a  neural  ' 

'  network  to  recognize  the  ship  silhouettes.  ' 

'  of  the  alphabet.  The  system  used  here  is  based  ’ 

'  on  a  two  layer  network  (not  including  the  input  ' 

’  layer)  with  20  neurons  in  the  hidden  layer  and  ' 

'  5  neurons  (one  for  each  ship  type)  in  the  output  ' 

'  layer.  The  moment  invariants  values  are  the  input  ' 

'  of  the  neural  network.  They  are  12  element  vectors' 

’  representing  the  invariants  for  the  silhouette.  ' 

'  The  network  has  already  been  trained  using  ' 

'  backpropagation  -  you  can  test  it  by  pressing  ' 

'  the  "New  Letter"  button.  This  passes  a  random  ’ 

'  letter  to  the  network.  The  "Noise"  slider  adds  ' 

'  random  noise  to  make  the  classification  problem  ' 

'  more  difficult.  The  "Rotation"  and  the  "Scale"  ' 

'  buttons  allovv  you  to  rotate  and  to  change  the  ' 

'  scale  of  the  silhouette  to  be  tested.  ' 

'  File  name:  interface. m 
helpwin (hlpStr , ttlStr) ; 


end; 


%  End  of  file  interface.m 
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%- 


Naval  Postgraduate  School  -  CA 


% 

Type  : 

Function 

% 

Name 

plotSilhouette .m 

% 

Function 

:  Draws  Ship  Sii' 

houetts  inside  t 

he  Graphical 

% 

Date  Cl 

march  20C1 

Version 

:  1.0 

% 

Author  ; 

Jorge  Amaral  Alves, 

LCDF  (Brazilian 

Navy) 

function  plotSilhouette (silhouette, rows, columns) 

Xl=zeros (rows , columns) ; 

XI ( : ) =silhouette; 

imagesc(Xl) ;%body  silhouette 

colormap (gray (2) ) 

set (gca, ’XTick‘ , [] , ’ YTick’ , [] ) ; 


%  End  of  file  plotSilhouette . m 
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%  Naval  Postgraduate  School  - 

%  Type  :  Function 

%  Name  :  segmentation .m 

%  Function  :  Segments  a  FLIR  real  image  using  a  histogram  and 
technique 

%  Date  01  march  2001 

%  Version  :  1.0 

%  Author  :  Jorge  Amaral  Alves,  LCDR  (Brazilian  Navy) 

%  =  =  =  ^  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

%Aircraft  Carrier 

[X,Map] =imread( ’ carrierlO -tif , ’ tif ’ ) ;%  240  lines  by  320  columns 
% [X,Mao] =imrsad ( ' carrierli . tif ‘ ' tif ' ) ; %  240  lines  by  320  columns 
%  [X, Map] =imread (' carriers . tif tif '); %  240  lines  by  320  columns 
%  [X,Map] =imread( ' carrierc . tif ,  ' tif ’ } ; %  240  lines  by  320  columns 

%De St rover 

% [X,Map] =imread( ’ab6 . tif ’ , ’ tif ’ ) ;%  240  lines  by  320  columns 
%  [X, Map] =imread ( ' ddl . tif ' ,  ' tif ' ) ; %  240  lines  by  320  columns 

%Msr chant 

%  [X, Map] =imread ( 'group3tanker.tif tif ■); %  240  lines  by  320  columns 
%  [X, Map] =imread( ’groups tankers .t if tif ’) ;%  240  lines  by  320  columns 
% [X, Map] =imread ( ' f lir2_2 -  jpg ' /  ' jp9 ' ) ; %  24  0  lines  by  32  0  columns 
%  [X, Map] =imread ( ' f lir2_3 . jpg ' ,  ' jpg ' ) ; %  240  lines  by  320  columns 
% [X,Map] =imread( ’ flir2_4 . jpg ’ , ’ jpg ’ ) ; %  240  lines  by  320  columns 
%  [X, Map] =imread ( ' f lir2_5 . jpg ' ,  ' jpg ’ } ; %  240  lines  by  320  columns 
% [X,Map] =imread { ' f lir2_6 . jpg ' , ' jpg ' ) ; %  240  lines  by  320  columns 
% [X,Map] =imread{ ’ f lir7_l . jpg ’ , ' jpg ’ ) ; %  240  lines  by  320  columns 
% [X, Map] ^imread ( ' f lir7_ll . jpg ’ , ' jpg ' ) ; %  240  lines  by  320  columns 
%  [X,Map] =imread ( ' f lir7_12 . jpg ‘ ' jpg ' > ; %  240  lines  by  320  columns 
% [X, Map] =imread ( ’ f lir8_l . jpg ' ,  *  jpg ' )  ; %  240  lines  by  320  columns 
%  [Z,Map] =imread { ' f lir8_12 . jpg ’ , ’ jpg ’ )  ; %  240  lines  by  320  columns 
% [X,Mapj =imread( ' f lir3_15 . jpg ' , ‘ jpg ' ) ; %  240  lines  by  320  columns 
%  [X, Map] =imread ( ' f lir3_2 -  jpg ’ ,  ' jpg ' ) ; %  24  0  lines  by  320  columns 

%  [Z, Map] =imread ( ' f lir7_5 . jpg ’ ,  ’ jpg ' ) ; %  24  0  lines  by  320  columns 

%Re search  Ship 

% [X,Mapj ^imreadC ' iranpbi . jpg ’ , ' jpg ' ) ; %  240  lines  by  320  columns 
%  [X, Map] “imread ( ' iranpb2 . jpg ’ ,  ' jpg ’ )  ; %  240  lines  by  320  columns 

%  [X,Map] =imread ( ' iranpbO . jpg ' ,  ' jpg ' ) ; %  240  lines  by  320  columns 

% fX, Map] =imread ( 'pointsur.bmp 'bmp '); %  240  lines  by  320  columns 


X1=X(:, :,1) ; 
figure 

colormap (gray (256)  ) 
imshow (XI ) 

[counts, X] =imhist (XI) ; 
figure 

stem (X, counts) 
p=polyf it (X, counts, 30) ; 
r=roots (p) ; 
maximo=max( counts) ; 
thresholds  f  ind  ( count  s ==inaximo ) 
rr= (r<threshold) .*r; 
newThreshold=abs (max (real (rr) ) )+10 
[x,y] =size (XI) ; 


threshold 
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XXl=:double  (XI)  ; 

Y=(XXl<newThreshold) .*(XX1>10) ; 

Y=double(Y) ; 

Y=Y*255; 

figure 

colormap (gray (256) ) 
imshow (Y) 

reslmg2  =  bwtnorph  (Y,  '  spur '  )  ; 
figure 

imshow  (res Img2)  ; 
colormap (gray (256) ) ; 

reslmg3  =  bwmorph (reslmg2 clean ') ; 
figure 

imshow(resImg3) ; 
colormap (gray (256) ) ; 

reslmg4  =  bwmorph (reslmg3 ,  ' f ill '  )  ; 
figure 

imshow (res Img 4 ) ; 
colormap (gray (256) ) ; 

[resImgS  ^  mapsl] =bwlabel (reslmg4 , 4 )  ; 
resImg6=zeros (size (Y) ) ; 
maximo=0 ; 
for  i=l:mapsl 

tempMap= (reslmg5==i) ; 
total=sum(sum(tempMap) ) ; 
if  ( total >maximo) 
resImg6=tempMap; 
maximo= total ; 

end 

end 

figure 

colormap (gray (256) ) 
imshow (res Img 6) 

resimg?  =  bwf ill (reslmg6 holes ') ; 
figure 

imshow (resimg?) ; 
colormap (gray (256) ) ; 

%Original  is  240x320  but  I  need  to  save  210x290 
Xl=reslmg7(16: (x-15) ,41:y) ; 

[x,y] =si2e (XI) ; 

abl=zeros (x*y, 1) ; 
abl=reslmg7 ( : ) ; 

save  abl  abl; 

%  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  :=  =  =  =  =  =  =  =  =  =  =  ==  =  =  =  =  =  =  =  =  =  = 
%  End  of  file  segmentation .m 
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Postgraduate  School 


%  Type  :  Main  program 

%  Name  :  createTestSet  .m 

%  Function  :  Creates  05  mat  files  containing  the  silhouettes 

%  of  each  ship  for  increments  of  one  degree  in  azimuth  and 

elevation 

%  then  plots  the  errors 

%  Date  01  march  2001 

%  Version  :  1.0 

%  Author  :  Jorge  Amaral  A^lves,  LCDR  (Brazilian  Navy; 


azAngleInc=l;%5 ;  %  ==>  36  azimuth  angles 
eleAngleInc=15 ; %  ==>  10  elevation  angles 
totalAzimuth=:180/azAngleInc;  %  =  =  >130 
totalElevation=2%4 ; %10 ; 

totalColumns=totalAzimuth*totalElevation; %each  column  means  one  view  ==> 
36*10=360  silhouettes 

totalLines=12 ;  %number  of  Tnoment  functions 

rows=210 ; 
columns=280 ; 

resultSet=zeros (totalAzimuth, totalElevation) ; 

1 

%aircraft  carrier 
[verts, faces] =aircarrier; 
rightResult=l ; 

resultSet=  f indResultSet (  verts, faces, azAngleInc,eleAngleInc,rightResult) ; 

aircarrierResultSet=resultSet ; 

save  aircarrierResultSet  aircarrierResultSet ; 

clear  aircarrierResultSet; 

2 

% destroyer 

[verts, faces] =destroyer ; 
rightResult=2 ; 

resultSet=  f indResultSet (  verts, faces, azAngleInc,eleAngleInc,rightResult) ; 

destroyerResultSet=resultSet ; 

save  destroyerResultSet  destroyerResultSet ; 

clear  destroyerResultSet; 

3 

% frigate 

[verts, faces] =f rigate; 
rightResult=3 ; 

resultSet=  f indResultSet (  verts, faces, azAngleInc,eleAngleInc,rightResult) ; 

f rigateResultSet=resultSet ; 

save  f rigateResultSet  f rigateResultSet ; 

clear  f rigateResultSet ; 

4 

%point  sur 

[verts , faces] =pointsur ; 


94 


rightResult=4 ; 

resultSet=  f indResultSet (  verts, faces, azAngleInc , eleAngleInc , rightResult) 

pointsurResultSet=resultSet ; 

save  pointsurResultSet  pointsurResultSet ; 

clear  pointsurResultSet; 

5 

%merchant 

[verts, faces] =merchant ; 
rightResult=5 ; 

resultSet=  f indResultSet (  verts, faces , azAngleInc , eleAngleInc, rightResult ) 

merchantResultSet=resultSet ; 

save  merchantResultSet  merchantResultSet ; 

clear  merchantResultSet; 

pause 

load  aircarrierResultSetl 
temp=aircarrierResultSet ; 
load  aircarrierResultSetl 
temp= [temp, aircarrierResultSetl ; 
ai rcarrierResult Set = temp; 

save  aircarrierResultSet  aircarrierResultSet ; 

load  destroyerResultSetl 
temp=destroyerResultSet ; 
load  destroyerResultSetl 
temp= [temp, destroyerResultSetl ; 
dest royerResult Set = temp; 

save  destroyerResultSet  destroyerResultSet ; 

load  frigateResultSetl 
temp=f rigateResultSet ; 
load  frigateResultSetl 
temps:  [temp,  frigateResultSetl  ; 
f rigateResul t Set = temp; 

save  f rigateResultSet  f rigateResultSet ; 

load  pointsurResultSetl 
temp=pointsurResultSet ; 
load  pointsurResultSetl 
temp= [temp, pointsurResultSetl ; 
pointsurResultSet=temp; 

save  pointsurResultSet  pointsurResultSet; 

load  merchantResultSetl 
temp=merchantResultSet ; 
load  merchantResultSetl 
temp= [temp, merchantResultSet] ; 
merchantResultSet=temp ; 

save  merchantResultSet  merchantResultSet ; 

%aircraft  carrier 
figure 

azimuthAngle= [-SO : 1 : 89]  ; 

errorAircarrier=sum (aircarrierResultSet ' ) ; 
errorPercent= (errorAircarrier/4 ) *100; 
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newAzimuthAngle=  [-90, ~85, -80, -75, -70, -65, -60, -55, ~50, -4  5, -30, - 
15,  0, 15,30,45,50,55,  60,65,  70,  75,80,  85]  ; 
newErrorPercentl=errorPercent (1:5:44) ; 
newErrorPercent2=errorPercent (45 : 15 : 136)  ; 
newErrorPercent3=errorPercent (140:5:179)  ; 

newErrorPercent= [newErrorPercentl , newErrorPercent2 , newErrorPercent3 ] ; 

plot (newAzimuthAngle, newErrorPercent , 'or'  ) 
hold 

plot  (azimuthAngle,  errorPercent) 
xlabel (' azimuth  angle') 
ylabel ( ' error (%) ' ) 

title { 'Aircraft  Carrier  Testing  Results  (adding  all  4  elevation  angles)') 
legend (' Training  set  ==>  96', 'Testing  set  -=>  720 Silhouettes  correct 

classified  ==>  642  (89.2%)')/ 

figure 

elevationAngle=  [0 : 15 : 46] ; 
errorAircarrier=sum(aircarrierResultSet) ; 
errorPercent = (errorAircarrier/180) *100; 
newElevationAngle= [0 : 15 : 46] / 
newErrorPercent=errorPercent (1:1:4) ; 
plot (newElevationAngle , newErrorPercent , ' or ' ) 
hold 

plot ( elevat ionAngle , errorPercent ) 
xlabel (' elevation  angle') 
ylabel ( ' error ( % )  '  ) 

title {' Aircraft  Carrier  Testing  Results  (adding  all  180  azimuth  angles)') 
legend (' Training  set  ==>  96', 'Testing  set  ==>  720 Silhouettes  correct 

classified  ==>  642'); 


%Destroyer 

figure 

azimuthAngle= [-90:1:89] ; 
errorDestroyer=suTn  (destroyerResultSet ' )  ,- 
errorPercent= (errorDestroyer/4 ) *100; 


newAzimuthAngle= [-90,-85,-80,-75,-70,-65,-60,-55,-50,-45,-30,- 
15,0,15,30,45,50,55,60,65,70,75,80,85] ; 
newErrorPercentl=errorPercent (1:5:44) ; 
newErrorPercent2=errorPercent (45:15:136) ; 
newErrorPercent 3 =errorPercent (140:5:179) ; 

newErrorPercent^ [newErrorPercentl , newErrorPercent2 , newErrorPercent 3]  ; 

plot (newAzimuthAngle, newErrorPercent , 'or' ) 
hold 

plot  (azimuthAngle,  errorPercent) 
xlabel ( ' azimuth  angle ' ) 
ylabel ( ' error ( % ) ' ) 

title (' Destroyer  Testing  Results  (adding  all  4  elevation  angles)') 

legend (' Training  set  ==>  96', 'Testing  set  ==>  720 ',' Silhouettes  correct 

classified  ==>  698  (96.9%)'); 

figure 

elevat ionAngle=  [0 : 15 ; 46] ; 
errorDestroyer=sum (destroyerResultSet) ; 
errorPercent = (errorDestroyer/180) *100; 
newElevationAngle= [0 : 15 :46]  ; 
newErrorPercent=errorPercent (1:1:4) ; 
plot (newElevationAngle, newErrorPercent , 'or' ) 
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hold 

plot ( elevationAngle , errorPercent ) 
xlabel {' elevation  angle') 
ylabel ( ' error { % } ’ ) 

title (' Destroyer  Testing  Results  {adding  all  180  azimuth  angles) ') 
legend (' Training  set  ==>  96', 'Testing  set  ==>  720 Silhouettes 

classified  ==>  698'); 


%?rigate 

figure 

a2imuthAngle= [-90 : 1 : 89]  ; 
errorFrigate=sum (f rigateResultSet ' ) ; 
errorPercent= (errorFrigate/ 4 ) *100 ; 

newAzimuthAngle= [-90,-85,-80,-75,-70,-65,-60,-55,-50,-45,-30,- 
15,0,15,30,45,50,55,60,65,70,75,80,85] ; 
newErrorPercentl=errorPercent (1:5:44) ; 
newErrorPercent2=errorPercent (45:15:136) ; 
newErrorPercent3=errorPercent  ( 140 :  5  : 179)  ,- 

newErrorPercent= [newErrorPercentl , newErrorPercent2 , newErrorPercent3 ] ; 

plot (newA2imuthAngle,newErrorPercent , 'or' ) 

hold 

plot (azimuthAngle, errorPercent) 
xlabel (' azimuth  angle') 
ylabel ( ' error { % ) ’ ) 

title (' Frigate  Testing  Results  (adding  all  4  elevation  angles)') 
legend (' Training  set  ==>  96', 'Testing  set  ==>  720 Silhouettes 

classified  ==>  666  (92,5%)'); 

figure 

elevationAngle= [0:15:46] ; 
errorFrigate=sum (f rigateResultSet ) ; 
errorPercent= (errorFrigate/180) *100 ; 
newElevationAngle= [0:15:46] ; 
newErrorPercent=errorPercent (1:1:4) ; 
plot (newElevationAngle, newErrorPercent , 'or' ) 
hold 

plot (elevationAngle , errorPercent) 
xlabel (' elevation  angle') 
ylabel ( ' error (%) ' ) 

title (' Frigate  Testing  Results  (adding  all  180  azimuth  angles) ') 
legend (' Training  set  ==>  96', 'Testing  set  ==>  720 Silhouettes 

classified  ==>  666'),- 


%Pointsur 

figure 

azimuthAngle= [-90:1:89] ; 
errorPointsur=sum (pointsurResultSet ' ) ; 
errorPercent=:  (errorPointsur/4 )  *100; 

newAzimuthAngle= [-90,-85,-80,-75,-70,-65,-60,-55,-50,-45,-30,- 
15,0,15,30,45,50,55,60,65,70,75,80,85]  ; 
newErrorPercentl=errorPercent (1:5:44) ; 
newErrorPercent2=errorPercent (45 : 15 : 136) ; 
newErrorPercent 3 =errorPercent (140:5: 179) ; 

newErrorPercent= [newErrorPercentl , newErrorPercent2 , newError Percents] ; 

plot (newAzimuthAngle, newErrorPercent, ' or ' ) 

hold 

plot  (azimuthAngle,  errorPercent) 


correct 


correct 


correct 
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xlabel ( ' azimuth  angle ' ) 
ylabel ( ’ error ( % )  ’ ) 

title (' Pointsur  Testing  Results  (adding  all  4  elevation  angles}') 
legend {' Training  set  =  =  >  96',  'Testing  set  ==>  720 Silhouettes 

classified  ==>  644  (89.4%)'); 

figure 

elevationAngle= [0:15:46] ; 
errorPointsur=sum (pointsurResultSet ) ; 
errorPercent= (errorPointsur/180 ) *100; 
newElevationAngle= [0 : 15 :46] ; 
newErrorPercent=errorPercent (1:1:4) ; 
plot (newElevationAngle , newErrorPercent , ' or ' ) 
hold 

plot (elevationAngle, errorPercent) 
xlabel (' elevation  angle') 
ylabel ( ' error ( % ) ' ) 

title (' Pointsur  Testing  Results  (adding  all  180  azimuth  angles)') 
legend (' Training  set  ==>  96', 'Testing  set  ==>  720 Silhouettes 

classified  ==>  644'); 

%M0r chant 
figure 

azimuthAngle= [-90 : 1 : 89]  ; 
errorMerchant=sum{merchantResultSet ' ) ; 
errorPercent= (errorMer chant /4) *100; 

newAzimuthAngle= [-90,-85,-80,-75,-70,-65,-60,-55,-50,-45,-30,- 
15,0,15,30,45,50,55,60,65,70,75,80,85]; 
newErrorPercent l=errorPercent  (1:5:44)  ,- 
newErrorPercent2=errorPercent (45 : 15 : 136) ; 
newErrorPercent3=errorPercent (140:5: 179) ; 

newErrorPercent= [newErrorPercentl , newErrorPercent2 , newErrorPercent3 ] ; 

plot (newAzimuthAngle , newErrorPercent , ’ or ' ) 

hold 

plot (azimuthAngle, errorPercent) 
xlabel (' azimuth  angle') 
ylabel ( ' error ( % ) ' ) 

title (' Merchant  Testing  Results  (adding  all  4  elevation  angles)') 
legend (’ Training  set  ==>  96', 'Testing  set  ==>  720 ',' Silhouettes 

classified  ==>  634  (88.1%)'); 

figure 

elevationAngle=  [0:15:46]  ; 
errorMerchant=sum (merchantResultSet ) ; 
errorPercent = (errorMerchant/180) *100; 
newElevationAngle= [0 : 15 : 46] ; 
newErrorPercent=errorPercent  (1:1:4)  ,- 
plot (newElevationAngle, newErrorPercent , 'or' ) 
hold 

plot (elevationAngle, errorPercent) 
xlabel (' elevation  angle') 
ylabel ( ' error ( % ) ’ ) 

title (' Merchant  Testing  Results  (adding  all  180  azimuth  angles)') 
legend (' Training  set  ==>  96', 'Testing  set  ==>  720 ',' Silhouettes 

classified  ==>  634'); 

%  End  of  file  createTestSet .m 


correct 


correct 


correct 


correct 
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%  Naval  Postgraduate  School  -  CA 

%  Type  :  Procedure 

%  Name  :  f indResultSet  .m 

%  Functior  :  returns  a  vector  with  the  size  of  all  the  viewangles  being 
tested 

%  where  "1"  will  mean  misclassif ied  and  ”0"  will  mean  correct 

classified 

%  Date  01  march  2001 

%  Version  :  1.0 

%  Author  :  Jorge  Amaral  Alves,  LCDR  (Brazilian  Navy) 


function  resultSet  =  f indResultSet ( 

verts, faces , azAngleInc, eleAngleInc, rightResult) ; 

totalAzimuth=180/azAngleInc ; 
totalElevation=2 ; %4 ; %10 ; %90/eleAngleInc ; 

totalColumns=totalAzimuth*totalElevation; %eac.h  column  means  one  view 
totalLines=12 ;  ^number  of  moment  functions 

resultSet=  zeros (totalAzimuth, totalElevation) ; 

trainingSet=  zeros (totalLines , 1) ; 

rows=210 ; 
columns=280 ; 

silhouettes=zeros  (rows*  columns,  total  Columns)  ; 


figNumber= figure (  ... 

' Name ' , ' Silhouette ' ,  ... 

'  Position *,  [120  120  280  210]);%  matriz  =>  210rov.'s  X  280columns 
f igure ( f igNumber ) 

patch ( 'Vertices ' , verts, ’ Faces ' , faces) ; 

view (3 ) 

axis  equal; 

axis  off 

axis  vis3d 

lineCount=l ; 

load  net; 

azimuthCount=0 ; 

for  i=-90 :azAngleInc : 89  %azimuth  popa  until  proa 

i 

azimuthCount=:azimuthCount+l  ; 
elevationCount=0 ; 

%for  j =0 : eleAngleInc : 16%46  %elsvation 

for  j =30 : eleAngleInc : 46  %elevation 

j 

e 1 e va t i onCoun t = e 1 e va  t i onCoun t + 1 ; 
view(i, j ) ; 

[X,map] =capture;  %[65,  65,  ... 

Xl=X-65;  %[0  0  ...  ]  max=l 

trainingSet (1 : 6 , lineCount) =f ind__mom_functions (Xl) ;  %  solid  silhouette 

X2=X1*255;  %[  0  0  ...]  max=255 
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fX'> 


XXl=edge (X2, 'prewitt ' ) ;  %only  the  edges  [0  0 

0  .  .  .  ]  inax==  1 

trainingSet (7:12, lineCount) =find_mom_f unctions (XXI) ;  %boundary 

A  =  sim(net, trainingSet) ; 
output=compet (A) ; 
result=f ind (output==l) ;%-l; 
result 

if  result==rightResult 

resultSet (azimuthCount , elevationCount)  =0; 
else 

resultSet (azimuthCount , elevationCount) =1 ; 

end 

end 

end 


End  of  file  f indRe suit Set .m 
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